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This Thesis is devoted to the study of the phenomenology of the Higgs sector of the 
minimal B — L extension of the Standard Model at present and future colliders. Firstly, 
the motivations that call for the minimal B — L extension are summarised. In addition, 
the model is analysed in its salient parts. Moreover, a detailed review of the phenomeno- 
logical allowed Higgs sector parameter space is given. Finally, a complete survey of the 
distinctive Higgs boson signatures of the model at both the Large Hadron Collider and 
the future linear colliders is presented. 
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Introduction 



1.1 The Standard Model 



Currently, the Standard Model (SM) of the electroweak (EW) and strong interactions 
of elementary particles represents one of the highest achievement of human kind in 
understanding the fundamental laws of Nature. 

The EW theory is based on the gauge symmetry group SU(2)l x U(1)y of lefthanded 
isospin and hypercharge. The Quantum Chromo-Dynamics (QCD) 1 , instead, is based 
on the symmetry group SU(3)c- These two theories provide a beautifully consistent 
picture of the particle physics phenomena observed up to now. 

The gauge symmetry groups of the SM give raise to a quantum field theory that is 
perturbative at sufficiently high energies and renormalisable. For a complete review of 
the model in its deepest theoretical and phenomenological implications, see [1, 2]. 

Starting from this framework, it is widely accepted that each SM particle (in particular, 
each massive gauge boson) obtains mass by mean of the so-called Higgs mechanism, i.e., 
the mechanism of spontaneous EW symmetry breaking (EWSB) realised by adding a 
theoretically consistent Higgs field: the minimal 2 choice is represented by the introduc- 
tion of a SU(2) doublet of complex scalar fields (see [3]). 

After the EWSB, the SU(2)l x U(1)y symmetry is spontaneously broken to the elec- 
tromagnetic U(1)q symmetry (where Q is the electric charge quantum number). Three 
of the four degrees of freedom of the doublet scalar field turn out to be absorbed in the 
longitudinal polarisation component of each of the three weak gauge bosons, W^ and 
Z, whilst the fourth one is the physical SM Higgs state h. At this stage, the fermion 

lr The theory of the strong interactions between the colored quarks. 
2 It respects the requirement of renormalisability. 
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masses are generated through the Yukawa interaction with the same scalar field and its 
conjugate one. 

Although the phenomenological power of this theory is remarkable (it has been deeply 
tested in various experimental scenarios), there is still no (either direct or indirect) 
evidence of the scalar sector impact in the probed phenomenology, in other words the 
Higgs boson has not been observed yet. 

The urgency of discovering the Higgs boson (and consequently establishing the nature 
of EWSB) has driven the Physics community to devote major efforts to the realisation 
of more and more powerful accelerators. 

For this, since decades a large number of machines (as LEP, SLC, Tevatron, etc.) has 
tried to investigate the SM to the high-precision measurement level (per mille accuracy), 
and apart from the astonishing agreement with the theoretical prediction concerning the 
fermion and vector boson sectors of the theory, the Higgs boson (the most important 
ingredient of the theoretical picture) is currently missing. 

At the same time, the recent observation of the pattern of neutrino oscillations (for a 
complete review, see [4-6]) produced the evidence of the inadequateness of the SM (in 
its minimal version) in describing the mass properties of neutrinos. 

Besides, dark matter {DM) evidence and related cosmological observations (see [7]) gave 
an hard blow to the minimal SM that, as well as for the massive neutrino case, can not 
be considered as a satisfactory theoretical framework for these phenomena. 

Other than by experimental facts, the SM is also affected by several theoretical prob- 
lems: firstly, quantum gravity is manifestely not included; secondly, the theory is affected 
by the so-called "Hierarchy problem" (see [8, 9]). 

However, while finding a solution to the former could be "postponed" to energy scales 
that are close to the Plank-scale O(10 19 )GeV , the latter is an indication of consistency 
problems already arising at the TeV energy scale. 

All these aforementioned theoretical and experimental reasons call for some extension 
of the SM. 

In fact, one could procede by extending the minimal SM by either top-down or bottom- 
up approaches. The former consist in the formulation of a theory that solves all the 
SM issues and can be broken down to a symmetry pattern that finally includes the 
SM as effective low-energy theory: a typical example of this approach is represented by 
supersymmetric scenarios as the (Next to) Minimal Supersymmetryc Standard Model 
(N)MSSM, see [10] ([11]). 
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The latter approach consists in piecing together a system of minimal extensions to give 
rise to the grandest possible system. In the context of this work, we choose to follow 
this approach in realising the so-called minimal B — L extension of the SM. 

1.2 A bottom-up approach: the minimal B — L model 

In order to realise a consistent extension of the SM, we exploit the fact that both the 
baryon (B) and the lepton (L) number are conserved quantities of the theory, as well as 
the B — L one. Furthermore, it is important to notice that the B — L number could be 
gauged in a U(l) symmetry group in combination with an augmented neutrino sector, 
and this creates a model that is free from anomalies. 

Hence, the minimal B — L extension of the SM consists of augmenting the SM gauge 
groups SU(3)c x SU (2)lxU (l)y by the aforementioned U(l) b-l symmetry (see [12, 13] 
and more recently [14-25]). 

This choice is minimal on three respects: 

• it is minimal in the gauge sector, enlarging the gauge group by adding one spon- 
taneously broken U(l) factor, which provides one new (neutral) gauge boson; 

• it is minimal in the fermion sector, adding one new heavy neutrino per generation; 

• it is minimal in the scalar sector, adding one new complex Higgs field, singlet under 
the SM gauge group. 

This extension gives rise to a model: 

• that is anomaly- free and gauge-invariant (see Section 2.1 for details); 

• that provides an "elegant" way to generate the light neutrino masses by means of 
the so-called "see-saw" mechanism (see [26-31] and Section 2.4 for details). 

This rather simple approach fulfils the phenomenological requirement of having a renor- 
malisable theory that provides a mechanism for giving mass to light neutrinos as well 
as a good candidate for DM (see [32, 33]). 

In addiction, it is important to notice that B — L symmetry breaking takes place at the 
TeV energy scale, hence leaving open the possibility of being part of a Grand Unified 
Theory (GUT) and giving rise to new and interesting TeV scale phenomenology. 
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In this work, we focus on the phenomenology arising from the main ingredients of EW 
plus B — L spontaneous symmetry breaking, i.e., the Higgs sector phenomenology of the 
minimal B — L model at present and future colliders. 

1.3 Experimental frameworks: present and future collid- 
ers 

In this Section we introduce the present and future machines that could possibly allow 
us to test the minimal B — L model. 

For this, two experimental frameworks have been taken into account: 

• the Large Hadron Collider (LHC). It is a hadron-hadron collider. It is the only 
high-energy accelerator that is currently working. Up to now, it is the world's 
largest and highest-energy particle accelerator ever made. 

• The International Linear Collider (ILC) and the Compact Linear Collider (CLIC). 
They are planned to be e + e~ colliders. They are two proposed accelerators and 
the projects have not yet been approved. If realised, they will be the largest and 
highest-energy linear particle accelerators ever made. 

1.3.1 The LHC 

The LHC was built by the European Organization for Nuclear Research (CERN) with 
the intention of testing various predictions of high-energy physics, with emphasis on 
addressing the nature of the EWSB. 

It is a circular collider and it is designed to investigate processes in which the initial 
state is characterised by either protons or heavy ions. 

It consists of several detectors/experiments: 

• ATLAS (A Toroidal LHC Apparatus): together with CMS is one of the two general 
purpose detectors; 

• CMS (Compact Muon Solenoid): together with ATLAS is one of the two general 
purpose detectors; 

• ALICE (A Large Ion Collider Experiment): it is designed to investigate the so- 
called "quark-gluon plasma" and the (de)confinement scenarios; 
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• LHCb (Large Hadron Collider b): it is mainly designed to investigate the nature 
of CP-violation in interactions of 6-hadrons. 

In this work we focus on discovery physics, i.e., we are interested on proton-proton 
collisions and possible minimal B — L signatures at both ATLAS and CMS. 

1.3.2 Future linear collider prototypes 

Although the LHC has finally entered its operational stage, a considerable part of the 
international physics community is focusing efforts in order to plan what the future of 
high-energy particle phenomenology and accelerators could be. 

Nowadays, the community is working on the proposal of two linear accelerator proto- 
types: the ILC and the CLIC. Both of them represent the new generation of electron- 
positron colliders. 

While the LHC and its multi-purpose detectors have been built with the aim of discov- 
ering new particles, the main goal that the new generation of linear colliders (LCs) is 
supposed to achieve is the subsequent profiling of the new physics whose evidence could 
arise in the next years at the LHC. 

Specifically, there are many phenomenological aspects on which the LHC is unsufficient 
in with respect to either ILC or CLIC: 

• if existing, direct measuring of mass, spin and coupling of the Higgs and new gauge 
bosons; 

• if existing, a complete profiling of any TeV scale extra-dimension; 

• if existing, a better understanding of any supersymmetric scenario and related 
dark matter aspects. 

In particular, in this work we are only interested in the minimal B — L model Higgs 
sector. Therefore we will consider the two linear accelerator configurations only for the 
scope of probing the existence of a Higgs boson. 

However, it is important to mention the fact that, if realised, the two accelerators will 
represent a unique way to precisely profile the Higgs sector of the minimal B — L model. 
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1.4 Organisation of the work 

This Thesis is organised as follows: 

• in Chapter 2 we describe the minimal B — L model in all its salient parts, with 
emphasis on the extended sectors, the B — L symmetry breaking mechanism, some 
formal aspects related to the gauge-invariance and ghost fields, and the "see-saw" 
mechanism for generating the light neutrino masses; 

• in Chapter 3 we perform a complete analysis of the experimental and theoretical 
constraints on the Higgs sector of the minimal B — L model, focusing on the bounds 
on the Higgs masses coming from direct searches at LEP, unitarity and triviality 
arguments; then, we also present the experimental and theoretical constraints on 
the other parameters of the model and discuss the "fine-tuning" in the minimal 
B — L model; 

• in Chapter 4 we present our phenomenological results of the investigation of the 
minimal B — L Higgs sector at colliders; we start with the study of Higgs branching 
ratios (BR), then we discuss some peculiar signatures at both the LHC and future 
LCs; 

• in Chapter 5 we summarise our result and we also conclude our work discussing 
some open issues of our research; 

• in Appendix A we list the potential of the minimal B—L scalar effective-Lagrangian; 
this part is useful in analysing formal aspects of theory, as it will be clear in Sec- 
tion 3.2; 

• in Appendix B we list the Feynman rules of the minimal B — L model adopting 
the Feynman-gauge. 



Chapter 2 

The minimal B — L model 



In this Chapter we introduce the minimal B — L model by a detailed description of the 
terms of its Lagrangian. 

As we have already intimated, this model is an extension of the SM, following the 
symmetry pattern: 

SU(3) C x SU(2) L x U(1) Y -> SU(3) C x SU(2) L x U(1) Y x U(1) B -l- (2.1) 

As we shall see, the charges of the two U(l) factors of the extended gauge group are 
associated with the SM weak hypercharge Y, and with the total (Baryon Lepton) 
number B — L 1 . 

The quantum consistency (anomaly cancellation) of the theory, together with the phe- 
nomenological requirement of generating neutrino masses, are fully satisfied by extending 
also the fermion content with a right-handed neutrino ur (singlet under the SM gauge 
group) for each family. 

Finally, the requirement of giving mass to the extra neutral gauge boson Z' is fulfilled 
by a minimal extension of the scalar sector as well, with the introduction of a scalar 
field x> singlet under the SM gauge group but not under the U{\)b-l (since it carries 
B — L charge). 

In Section 2.1 we begin by identifying the model Lagrangian in all its parts. In Section 2.2 
we then proceed with a discussion of the extended Higgs mechanism, breaking SU(2)l x 

In the general B — L model, this identification is not generally possible since the presence of "mixing" 
makes it basis-dependent, therefore the covariant derivative is generally non-diagonal. Indeed, it is 
always possible to diagonalize the covariant derivative by defining an "effective charge" Y p (with the 
corresponding coupling g p ), which is a linear combination of Y and B — L, but in this case a more 
careful treatment should be done and it lie outside of our purposes. 
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U(1)y x U(l) B—L down to U(1)q , and of the ensuing spectrum of gauge bosons, fermions 
and scalars. In Section 2.3 we present a brief analysis of the minimal B — L gauge-fixing. 
In Section 2.4 we introduce the principles of the so-called "see-saw" mechanism. Finally, 
in Section 2.5 we conclude this Chapter with a summary of the introduced parameters 
in comparison with the SM case. 

It is important to mention the fact that the content of this Chapter will allow one to 
get all the necessary elements to implement the minimal B — L model in any Feynman- 
rules generator (as LanHEP [34], FeynRules [35], etc.). In Appendix B we collect the 
Feynman rules of the model that we have obtained by mean of the LanHEP package. 

2.1 The Lagrangian and the parameterisation 

The minimal U(1)b-l extension of the SM is realised by augmenting the symmetry 
structure by a U(l) gauge group, that is related to the B — L symmetry. 

Without loss of generality, the classical gauge-invariant Lagrangian for the class of mod- 
els under consideration (obeying the SU(3)c x SU(2)l x U(l)y x U(1)b-l) can be 
decomposed as: 

if = S? 8 + 3?YM + S?f + £Y, (2-2) 

where the various terms represents the scalar, Yang-Mills, fermion and Yukawa part, 
respectively. 

The general structure is similar to the SM Lagrangian, although (due to the different 
symmetry structure) each term takes into account the differences in the gauge, fermion 
and scalar sectors. 

In the following, sector by sector, we will highlight analogies and differences with respect 
to the SM case. 

2.1.1 The Yang-Mills sector 

As in the SM, the vector fields are uniquely determined by the choice of the gauge group, 
and by the transformation in their adjoint representation. Hence, in the J^ym, the non- 
Abelian field strengths therein are the same (as in the SM) and the only difference is 
contained in the Abelian terms. 
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Explicitely, it is formalised as follows: 

&ym = -\G^ a CT a - \w^W^ a - \f^F, v - \f^F^ , (2.3) 

where 

F^ = d^Bv - d v B^ (2.4) 

K» = ^K-dvB'^, (2.5) 

here, B and B' are the gauge fields associated with U(l)y and U(1)b-l, respectively. 
2.1.2 The fermion sector 

The fermion content of the model is the same of the minimal SM , except for the addition 
of a right-handed neutrino vr (singlet under the SM gauge group) for each of the three 
lepton families. 

As we already mentioned, this addition is essential both for anomaly cancellation and 
preserving gauge invariance. 

Following the field-basis notation introduced in Subsection 2.1.1, the covariant deriva- 
tives of the B — L model are defined as the usual SM non-Abelian part plus an Abelian 
part: 

D„ = d„ + ig s T a G« + igT a W fl a + i 9l YB^ + i(gY + ^Yb-l)^. (2.6) 

In the minimal version of the B — L, g = is assumed (i.e., no mixing between the two 
U(l) factors) 2 and the covariant derivative becomes: 

D„ = d„ + ig s T a G« + itfFW* + i 9l YB^ + %^Y B -l%. (2.7) 

Then, the fermionic Lagrangian is given by: 

3 

= ^ {nkLl^D^qkL + iukR'y f _ l D fJl u kR + idkR^D^dkR + 
fc=i 

+ilk~L7^l kL + imil^e kR + iVk~Rl^v k Rj , (2.8) 

where the charges of the fields are the usual SM ones, plus the B — L ones: 

2 It is important to highlight that this condition holds at the EW scale only: if we assume the running 
of g, then it will monotonically grow spoiling the "minimality" of the model (see Subsection 3.2.2 for 
details) . 
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• Yb-l = 1/3 for quarks, 

• Yb-l = — 1 for leptons. 

Assuming that the conjecture of "universality" is true, no distinction between genera- 
tions has been made. 

2.1.3 The scalar sector 

The model under study has an extended gauge sector, with an additional neutral gauge 
boson Z' with respect to the SM. In order to realise a consistent Higgs mechanism 
(giving mass not only to the SM weak bosons but also to the Z') it is necessary to 
enlarge the SM Higgs sector by means of a further complex Higgs singlet X - 

The B — L charges of the two scalar fields are set to be: 

• Yb-l = 0, 

• Y*_ L = +2; 

this choice is essential to preserve the gauge invariance of the model. 
The most general gauge- invariant and renormalisable scalar Lagrangian is: 

JSf s = (ZW) f D^H + (D» X ) j D„ X ~ V(H, X ) , (2.9) 
with the scalar potential given by 

V(H, X ) = m 2 H^H + ^\x\ 2 + (H^H \ X f ) ^ ) ( ^ p ) 

= m 2 rfH + | x I 2 +\i(H^H) 2 + \ 2 | X I 4 +X 3 H^H \ X \ 2 ■ (2.10) 

2.1.4 The Yukawa term 

To complete the description of the Lagrangian in equation (2.2), it is necessary to per- 
form a treatment of the Yukawa couplings. 

In addition to the Yukawa couplings of the minimal SM, we have two new types of 
Yukawa interactions involving right-handed neutrinos: 

= -y < j k qLjd Rk H - yj k qEjU Rk H - y e jk l Lj e Rk H 

-yjk l LjVR k H - yf k (yR))v Rk x + h.c, (2.11) 
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where H = ia 2 H* and i, j, k run from 1 to 3. 

It is important to notice that the Yukawa interactions can generate both Dirac mass 
terms and Majorana mass terms for right-handed neutrinos (both of them in the secon 
line of equation (2.11)). As we will see in Section 2.4 these are the essential ingredients 
of the "see-saw" mechanism for giving mass to light and heavy neutrinos. 

As aforementioned, it is now clear, from the terms involving the x scalar field, that 
Y% L = +2 is needed in order to ensure the gauge invariance. 



2.2 Spontaneous SU{2) L x U(1) Y x U{1) B -l breaking 

We generalise the SM discussion of spontaneous EW SB to the more complicated case 
represented by the potential of equation (2.10). 

To determine the condition for V(H, x) to be bounded from below, it is sufficient to 
study its behaviour for large field values, controlled by the matrix in the first line of 
the aforementioned equation. Requiring such a matrix to be positive-definite gives the 
conditions: 

4AiA 2 -A| > 0, (2.12) 
Ai,A 2 > 0. (2.13) 

If the above conditions are satisfied, the choice of parameters is consistent with a well- 
defined potential, hence we can proceed to the minimisation of V as a function of constant 
Vacuum Expectation Values (VEVs) for the two Higgs fields. 

Since the minimisation procedure is not affected by the choice of the gauge, it is not 
restrictive to define the two VEVs in the following way: 



M S ^ (2.14) 
with v and x real and non-negative. 

Then, the search for extrema of V is made by mean of the following differential set of 
equations: 

9V, , ( 2 , 2 A 2 \ 

— (v.x) = v • m \\v + —x = 

9V, , / 2x 2 a! 2 \ [ ' 

— (v,x)= x • f /rA 2 ar + y v ) = 
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The physically interesting solutions are the ones obtained for v, x > 0: 



-A 2 m 2 + 4f ^ 



X^X 



1^2 



A 3 



AiA 



1^2 



4 



(2.16) 
(2.17) 



Since the denominator in equations (2.16)-(2.17) is always positive (assuming that the 
potential is well-defined), it follows that the numerators are forced to be positive in order 
to guarantee a positive-definite non- vanishing solution for v and x. 

In order to identify the extrema, we need to evaluate the Hessian matrix: 

/ d 2 v d 2 v\ 



H(v, x) 



dv 2 dvdx 
d 2 V d 2 V 
\ dvdx dx 2 ) 



2\\v X^vx 

X3VX 2\2X 2 



(2.18) 



From this equation, it is straightforward to verify that the solutions are minima if and 
only if equations (2.12)-(2.13) are satisfied. 



2.2.1 The scalar mass spectrum 

To compute the scalar masses, one must expand the potential in equation (2.10) around 
the minima found in equations (2.16)-(2.17). 

Since the physical mass eigenvalues are gauge invariant, we define the Higgs fields fol- 
lowing the unitary-gauge prescription: 



H = 





h + v ) , X ^^7T- 
V V2 



After standard manipulations, the explicit expressions for the scalar mass eigenvalues 
are: 



M 2 hi = X lV 2 + \ 2 x 2 - v^Aiv 2 - A 2 x 2 ) 2 + {Xzxvf, (2.20) 
M 2 2 = X lV 2 + X 2 x 2 + v^Ai^ 2 - A 2 x 2 ) 2 + (\ 3 xv) 2 , (2.21) 

where hi and h 2 are the scalar fields of definite masses and M^ 2 respectively, and 
we conventionally choose < M% . 
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These eigenvalues are related to the following eigenvectors: 



hi 

h 2 



cos a — sin a \ I h 
sin a cos a I \ h' 



(2.22) 



where -f < a < § fulfils 3 : 



sin 2a 



cos 2a 



^(Xiv 2 - X 2 x 2 ) 2 + (A 3 xu) 2 ' 

A2X 2 — Ait; 2 
x /(Ai^-A 2 x 2 ) 2 + (A 3 ^) 2 ' 



(2.23) 
(2.24) 



For completeness, it is useful to write the isomorphic transformation between the two 
A1-A2-A3 and M^-Mfo-a spaces. 

From equations (2.20)-(2.21)-(2.23), it is straightforward to have: 



Ai 
A 2 
A3 



H — sin a = — -± cos a + 4- sin a 



2?; 2 2u 2 

2x 2 2x 2 
[Ml 2 - Ml) 



2v 2 



2v 2 



2 M L . 2 , M h 2 2 

H r— ^ cos a = - n sm a + _ „ cos a 



2x 2 



2x 2 



sin (2a). 



(2.25) 



2.2.2 The gauge mass spectrum 

To determine the gauge boson spectrum, we have to expand the scalar kinetic terms as 
in the SM case. From consistency arguments, we expect that there exists a massless 
gauge boson (the photon) and the other gauge bosons acquiring mass. 

Moreover, we expect that the SM charged vector boson spectrum is not affected by the 
extension (since it involves the non-Abelian part only). As for the Abelian vector bosons, 
using the unitary-gauge parameterisation, we write the kinetic terms in equation (2.9) 
as: 



-dPhdph + -{h + v?{Q\) gW a »a a + gi B» 



1 



l^hdfji + l -{h + v) 2 [g 2 \W/ - iW 2 "\ 2 + (gW/ - 9l B») 2 ] , (2.26) 



3 In all generality, the whole interval < a < 2ty is halved because an orthogonal transformation is 
invariant under a — > a + n. 
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and 



^tid^ti + ±(h' + x) 2 (g[2B"*) 2 . 



(2.27) 



In equation (2.26) we can easily identify the SM charged gauge bosons W^, with M\y = 
gv/2. 

As for the other fields, the vanishing mixing between the two U{1) factors (i.e., the 
assumption g = 0) allow us to immediately identify both the SM-like pieces and the 
new vector boson Z' = B'. 

Once we set the B^, W% and B tiL as field basis, the explicit expression for the squared 
mass matrix is: 



M' 



v 

T 



o 



9x ~99i 
991 g 2 o 

2 



V o o i6$K) 2 y 

( sin 2 $ w — cos $ w sin i? w 
cos i?,„ sin #,„ cos 2 d w 




-(9 2 +9i: 



\ 










16x 2 (g[) 2 

v 2 (g 2 +g()J 



where we have made use of the well-known relations: 



cos i?, 



VtFTgJ' 



sin 



91 



Vg^TgJ 



(2.28) 



(2.29) 



As in the SM case, if we want to diagonalise the 2x2 sub-matrix in equation (2.28), 
we need to apply a rotation along the B'(Z') field, defined by: 

RewW w 



^cos $ t 
sin 

V 

and this allow us to isolate each mass eigenvalue: 



- sin -d w o\ 
cos ~& w 
1/ 



Rew^M^RewWw)]' 1 







W + <?f) 
\0 4x 2 ( 5 ' 1 ) 2 y 



(2.30) 



(2.31) 
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Finally, we can associate the mass eigenvalues with the corresponding physical vector 
boson eigenstates: 



2.3 The Feynman-gauge and the ghost Lagrangian 

Once the Lagrangian is established and the gauge symmetries are spontaneously broken, 
the gauge bosons acquire their masses and we have constructed a consistent theory. 

However, the way we treated the Higgs mechanism so far is too simplicist: assuming 
the unitary-gauge framework, we have not depicted any exhaustive description of the 
computational Feynman rules of the minimal B — L model, hence we are not ready for 
any phenomenological analysis yet. 

Moreover, as we shall see in Subsection 3.2.1, we are missing a correct treatment of 
formal unitarity aspects of the theory. 

Finally, it is well known that if the definition of each model is delivered in the Feynman- 
gauge the public softwares dedicated to the computation of amplitudes/cross-sections 
(as CalcHEP [36, 37], MadGraph [38], FeynArts-FormCalc [39], etc.) increase their 
computational power, hence with an explicit parametrisation of the Goldstone contri- 
bution within the Higgs potential. 

For this, in the Subsection 2.3.1 we will introduce a standard parametrisation for the 
Higgs fields in the Feynman-gauge. 

As for the aforementioned unitarity of the theory, it will be apparently spoiled by the 
fact that the Feynman gauge will introduce the effects of unphysical particles in the 
computational formalism. In order to cancel this effects, in the Subsection 2.3.2 we will 
define a gauge-fixing Lagrangian, using the Fadeev-Popov method, in order to restore 
the computational consistency of the the minimal B — L model. 

2.3.1 The Feynman gauge 



M A = 




(2.32) 



M z , = 2xg' 1 



We focus again on the scalar sector of the Lagrangian described in Subsection 2.1.3. 



Chapter 2. The minimal B — L model 



16 



Following the prescription of the Feynman-gauge parametrisation of the Higgs doublet 
and singlet of the minimal B — L model, we consider not only the Higgs fields and their 
VEVs, but also the so-called Goldstone bosons. 

Hence, H and x are now defined in the following way: 



where = w 1 ^iw 2 , z and z' are the Goldstone bosons of , Z and Z\ respectively. 

From this definition and the set of relations described by equations (2.20)-(2.25), we can 
calculate the explicit rules of the interactions (i.e., the Feynman rules) in terms of mass 
eigenstates and couplings. 

As we will show in Subsection 3.2.1, a particular interesting aspect of the structure of 
the scalar Lagrangian in the Feynman gauge is that it fully describes the interaction 
of longitudinally polarised vector bosons and Higgs bosons in the high energy limit 
("equivalence theorem", see the Chapter 21 of [1] and [40]). 

Moreover, if we assume that the couplings of the theory are perturbative and small 
then we can also apply the prescription D„ ~ in order to obtain the would-be- 
Goldstone scalar sector: to all intents it works as an effective high energy theory of the 
longitudinally polarised vector bosons and Higgs bosons interactions. 

For this, the would-be-Goldstone effective theory is a particulary useful tool for analysing 
the perturbative unitarity stability of longitudinally polarised vector boson scatterings 
in the high energy limit, and we list the complete set of the functions appearing in the 
would-be-Goldstone Lagrangian, in Appendix A. 

2.3.2 The Fadeev-Popov Lagrangian 

Since there is no mixing in the neutral boson sector, the gauge-fixing in the minimal 
B — L model is fulfilled by a rather easy procedure. 

Firstly, the Higgs doublet gauge-fixing Lagrangian is not affected by the B — L extension. 
Hence, following the notation of [1], the ghost Lagrangian in the Feynman gauge is the 
usual iSM-like one: 




(2.33) 



= c 



(d^D fl ) ab -g 2 (T a (H))-T b H c b 



(2.34) 



where the c's are the gauge-fixing fields corresponding to the W's and Z gauge bosons. 



Chapter 2. The minimal B — L model 



17 



Due to the fact that the x held belongs to an Abelian gauge symmetry structure, we 
have to extract the gauge-fixing for the Abelian case in the Feynman gauge. 

To begin, we focus on the spontaneously broken term of the scalar Lagrangian in equa- 
tion (2.9): 

J? SAb = {D»x) ] D„ X ~ V(H, X ), (2.35) 
with = 8^ + 2ig' 1 Z'. From the second equation of (2.33), we obtain: 



&.a> = - Ig'iZ'^') 2 + \{d,z> + 2g> 1 Z^x + h')) 2 - V(H, X ). 



(2.36) 



This Lagrangian is invariant under an exact local symmetry: 



Sti = -a(x)z', 5z' = a(x)(x + ti), SZ 1 = -—jd^x). (2.37) 

2<7i 



Due to the existence of this local symmetry, in order to define the functional integral 
over the variables h' , z', Z', we must introduce the Fadeev-Popov gauge-fixing. 

Following the standard techniques (see the Chapter 9 of [1]), we define the functional 
integral: 

(2.38) 



Z = j V{Z')V(h')V{z') exp i J Sf t 



Introducing a gauge-fixing constraint we find: 
Z = C J V(Z')V(h')V(z') exp % J S?s A 



5(G{Z',h',z'))det 



6G 
5a 



(2.39) 



where C is a constant proportional to the "volume" of the gauge group and G is a 
gauge-fixing condition. At this point, we can introduce the gauge-fixing constraint as 
5(G(x) — u}(x)) and integrate over uj(x) with a Gaussian weight, to obtain: 



Z = C J V(Z')V{h')V(z') exp i J d 4 x(J? SAb - ^G 2 ) 



(2.40) 



The gauge-fixing function G is arbitrary, and a common choice in the Feynman-gauge 
is: 

G = d»Z' - 2g[xz'. (2.41) 



From the gauge-fixing condition, it is straightforward to obtain the Lagrangian of ghost 
for the Abelian terms of the scalar Lagrangian. Firstly, we evaluate the gauge variation 
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of G: 



5G 
5a 



2<?i 



T d 2 - 2g' lX (x + ti) 



1 



-d 2 - (2g[x) 2 1 + 



h! 



(2.42) 



The determinant of this operator is related to the ghost Lagrangian by the equation: 

'SG 



where 



det 



5a 



VcVc exp 



2.9 



,5G 

5a 



„z' _ -z> 



i / d xJfpp 



1 + ^1 



Z' 



(2.43) 



(2.44) 



In the last equation, Mz< = 2g[x and c 's are the ghost fields related to the Z' boson. 

It is interesting to notice that since this belongs to an Abelian gauge structure, the ghost 
fields does not couple directly to the gauge field, but only to the physical Higgs field. 



2.4 See-saw mechanism and neutrino masses 

In the SM framework, there is not a straightforward way to generate the experimentally 
observed neutrino masses and oscillations: any isolated extension that solve this open 
issue (effective Majorana mass terms, sterile right-handed neutrinos, etc.) is affected by 
consistency problems (they spoil either renormalisability or gauge invariance, see [4] for 
details). 

The minimal B — L model provides an elegant "natural" solution: the presence of right- 
handed neutrinos in the Yukawa Lagrangian (equation (2.11)) gives raise to the so-called 
"see-saw" mechanism. 

In details, after the spontaneous breaking of the gauge symmetry, the Dirac neutrinos 
combine to six Majorana mass eigenstates, and the mass matrix is: 



Ji=\ \ , (2.45) 




where mo and M are respectively the Dirac and Majorana mass matrices, defined as: 

m D = %Lu, M = V2y M x. (2.46) 
v2 

Once assumed that the hierarchy Ad <C Am (where A indicates the energy scale) is true, 
the diagonalization of the mass matrix realises the "see-saw" mechanism. 
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After this procedure, we are left with three light Majorana neutrinos v\ and three heavy 
Majorana neutrinos v^, whose 3x3 mass matrices, denoted by Mi and Mh respectively, 
are given by: 

M h ~ M = V2y M x. (2.47) 



In equations (2.47) we can appreciate the "see-saw" effect: the greater is M, the smaller 
is M t . 

As explicit example we consider the case in which there is no mixing between the three 
generations, hence the matrices mo and M are two diagonal matrices: in this case the 
6x6 matrix splits in three 2x2 matrices. 

Thereafter, the diagonalisation is realised by the transformation 

/ cos OLi — sinaA / m' n \ / cos a; sinaA / (mV>) 2 A 

-diag -i-^-,M l , (2.48) 



\sinaj cosaj J \m D M l J \ — sinaj cos on J V M 

where i = 1,2,3 denotes one of the three generations and ot = arcsin (m l D /M 1 ^. 

We can roughly estimate the mass-scale Am needed to obtain neutrino masses in agree- 
ment with phenomenological constraints [41, 42]. 

By taking A/ < 1 eV and Ad ~ EW scale we get: 

A 2 

Ai ~ — < 1 eV => A M > 10 13 GeV. (2.49) 
Am 



We should keep in mind, however, that A^> could well be several orders of magnitudes 
smaller than the weak scale (for example, the electron mass): in such a case, much 
smaller scales for Am are allowed. 

Anyway, starting from equations (2.46), a generalised condition could be: 

v\y v \ <x|y M |. (2.50) 



It is important to recall that in the general case, in order to extract the individual mass 
eigenvalues and eigenstates of the mixed neutrinos eigensystem, we must separately 
diagonalise the two mass matrices M; and M^. In the case of Mi, the diagonalisation 
gives rise to the well-known mixing Upmns ( see [43-46]) matrix with 6 real independent 
parameters, 3 angles and 3 phases. 
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2.5 Counting of parameters 

The SM contains 18 real parameters (3 in the gauge sector, 2 in the Higgs sector, 13 in 
the fermion sector). 

Once we add neutrino masses to the SM via the see-saw mechanism, and assume no 
mixing in the right-handed neutrino sector 4 , 12 more real parameters are introduced (6 
neutrino masses and the 6 parameters of Upmns), f° r a total of 30. 

The Yang-Mills sector contains 1 more parameter than in the SM model: the gauge 
coupling constant g[. 

Finally, The scalar potential depends on 7 real parameters, the two scalar masses (to 2 , 
fi 2 ), the three couplings (Ai, A2, A3) and the two VEVs (v, x): three of them are equiv- 
alent to the corresponding SM parameters, thus we have here 4 additional parameters 
with respect to the SM. 

In conclusion, the total number of parameters is 35. 



4 This is realised by neglecting the additional complex parameters of their Yukawa couplings. Let's 
remind that in this work we are mainly interested in colliders and discovery physics, hence allowing any 
mixing either in the heavy or light neutrino sector does not affect our conclusions. 



Chapter 3 

The Higgs sector parameter space 



As we have already seen in Section 2.5, the B — L extension of the SM introduces a new 
set of parameters. 

Apart from the Higgs boson mass, all the parameters involved in the definition of the 
SM Lagrangian have been set by experiments; when we call for the minimal B — L 
extension, it opens a new parameter space defined by 12 parameters (from the neutrino 
sector 1 ), 1 parameter (from the B — L gauge boson sector 2 ) and 4 parameters (from the 
Higgs bosons sector 3 ). 

In principle, this could allow to explore potentially infinite new phenomenological im- 
plications, nevertheless one has to take into account the experimental and theoretical 
constraints that affects each sector. 

As for the Higgs sector, in the past four decades an enormous effort has been made to 
improve both the theoretical and experimental techniques that allow us to understand 
what is the forbidden parameter space. 

We can split these procedures in two sets: "experimental" and "theoretical" . 
The ones belonging to the former give constraints: 

• by EW precision tests (Subsection 3.1.1), 

• by direct searches (Subsection 3.1.2); 

the ones belonging to the latter give constraints: 

lr The 6 neutrino masses and the 6 parameters of Upmsn- 
2 The g[ coupling. 
3 The 3 A's and the x. 
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• by perturbative unitarity arguments (Subsection 3.2.1-3.2.3), 

• by triviality and vacuum stability (Subsection 3.2.2-3.2.3), 

• by "naturalness" arguments (Subsection 3.2.4). 

In this Chapter we present the main results that we have obtained by the application of 
these techniques to the minimal B — L model and their major implications. 

3.1 Experimental constraints on the Higgs boson sector 

Despite the fact that the Higgs boson(s) has(have) not been discovered yet, we know by 
consistency that the Higgs sector should participate to the known phenomenology. 

Firstly, the Higgs bosons (as well as the other new particles of this model: heavy neu- 
trinos and Z') must contribute to the quantum corrections to the high-precision EW 
observables 4 , and this could impose boundaries on the free parameters. 

Moreover, there are constraints coming from direct searches of Higgs bosons at colliders, 
in particular from LEP. 

In this Section we will summarise the main results of these indirect and direct (respec- 
tively) constraints. 

3.1.1 High-precision data and constraints on the Higgs sector 

Even if the B — L model is a minimal extension of the SM, it shows a phenomenological 
richness that imposes an accurate approach to the EW precision data analysis. 

In principle, if one assumes that the symmetry group of new physics (NP) is still 
SU(3)c x SU(2)l x U(1)y, it is possible to parametrise the radiative corrections in 
such a way that the contributions from NP could be easily implemented and con- 
fronted with EW precision data (a popular example is the well-known Peskin-Takeuchi 
parametrisation in [47, 48]). 

However, it is important to mention the fact that the parametrisation is based on the 
fact that the Higgs bosons must be the only source of quantum corrections. 

Obviously, this is not the case of the minimal B — L model, because of its extra U (1) 
extension. 

4 In the last twenty years an enormous amount of EW precision data have been collected by the e + e~ 
colliders LEP and SLC. 
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Indeed, in this case we have not only a new gauge boson Z' that affects the choice of 
the precision parameters parametrisation, but also the presence of heavy neutrinos that 
are not weakly coupled with the boson sector: this spoils the first assumption. 

In the absence of new gauge bosons or heavy neutrinos, many successful attempts have 
been made to profile the constraints on the Higgs sector parameter space (in the multiple 
Higgs scenario, some interesting study has been performed in [49]), but unfortunately 
the B — L model falls outside these conclusions. 

Hence, in this work we have not considered any boundary condition coming from the EW 
precision analysis. However, we borrow a qualitative conclusion from the aforementioned 
multiple Higgs scenario paper: we assume that the Higgs bosons inversion limit a — > ir/2 
is not allowed. 

3.1.2 Direct searches and constraints on the Higgs sector 

The minimal SM Higgs boson has been searched in the LEP experiments, firstly at 
energies close to the Z boson peak (LEP1), then with center of mass energies up to 
200 GeV (LEP2). 

The main production modes that have been explored at LEP1 are the so-called Bjorken 
process (Z — > HZ — > Hff) and the Z — > H^ (through triangular loops of t's and Ws). 

In order to avoid the large hadronic background, Higgs boson have been unsuccessfully 
searched in the two aforementioned channels (see [50-53]), leading to a 95% Confidence 
Level (CL) limit of Mh > 65 GeV on the SM Higgs mass. 

In the LEP2 energy regime, the Higgs boson have been searched in its dominant pro- 
duction process: the Higgs-strahlung, Z — > HZ. 

LEP collaborations have explored several topologies, combining their results in the anal- 
ysis of the Z boson recoil energy that led to the well known exclusion limit of the SM 
Higgs boson mass (see [54]): 

M H > 114.4 GeV (3.1) 
at the 95% CL from the non-observation. 

In the B — L model (as well as in other extensions of the SM), this limit is shifted back 
in mass by the fact that Higgs bosons have non-standard coupling to the Z (it is reduced 
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by cos a) : 



<7SM(e H 



HZ) 



a B ^ L (e+e- -> H X Z) 



cos 2 a. 



(3.2) 



Indeed, the LEP collaborations have also released a detailed analysis of the 95% CL 
exclusion limits on the (lightest) Higgs boson mass as function of the reduced coupling 
£ 2 = cos 2 a(sin 2 a) for the hi(fi2), and it is shown in figure 3.1. 

In this work, we will consider the latter as the only experimental constraint on the B — L 
Higgs sector parameter space. 



3.1.3 Other experimental constraints 

In view of an extensive numerical analysis of the phenomenological implications at col- 
liders, one has to have a complete understanding of how other experimental bounds 
could affect the B — L parameter space. 

For this, we briefly present the experimental constraints on the gauge boson and neutrino 
sectors that we are going to consider in the rest of this work. 
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• Mz>, the new gauge boson mass. An indirect constraint on Mz> comes from anal- 
yses at LEP of precision EW data (see [55], based on the analysis of experimental 
data published in [56-60] ) 5 : 

M 7 i 

-f- > 7 TeV . (3.3) 
9i 

Further limits have been obtained at Tevatron [20, 62, 63]. Both have been taken 
into account in this work. 

• m Ul , the SM (or light) neutrino masses. We use the cosmological upper bound 
^ m Vl < 1 eV [64]. Ultimately, they have been taken to be m Vl = 10~ 2 eV. Since 
we are mainly interested in the discovery and collider physics phenomenological 
aspects, we remind that in the rest of our work, for our illustrative purposes, we 
take all neutrino masses (both light and heavy) to be degenerate. 



3.2 Theoretical constraints on the Higgs boson sector 

Several theoretical methods have been developed to establish what are the constraints 
on the Higgs parameter space. 

In this Section we apply them to the minimal B — L model. 

They tipically follow from general principles, like the assumption that the perturbative 
unitarity and the couplings of a theory are preserved to some energy scale (Subsec- 
tion 3.2.1), as well as the stability of the symmetry breaking vacuum expectation value 
(Subsection 3.2.2). 

These arguments could also be extended to different sectors of the theory in order to 
constrain other free parameters, as g[, as we will show in the Subsection 3.2.3. 

Finally, one could also propose some conjecture about the "naturalness" of the model, 
finding reasonable to indicate some "natural" range of value for the allowed parameter 
space of the theory in order to minimise the so-called "fine-tuning" problem (Subsec- 
tion 3.2.4). 

3.2.1 Perturbative unitarity 

Even if there is no evidence of the (B)SM Higgs boson(s), the need of such kind of 
particle in order to guarantee the perturbative unitarity of a theory is of fundamental 
concern. 

A less conservative approach, based on Fermi-type effective four-fermions interactions, gives the 
weaker constraint — f- > 6 TeV [61]. 
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In fact, the potential problem that affects the SM and its extensions is the following: in 
the absence of the Higgs boson contributions, the longitudinally polarised vector boson 
interactions violate unitarity at some energy scale. 

To see this, we must recall some general properties of the scattering amplitudes; firstly, 
the amplitude A of a 2 — > 2 process could be decomposed into partial waves ai of orbital 
angular momentum I: 

oo 

A = 16irJ2(2l + l)P l (cosQ)ai, (3.4) 

1=0 

where Pi are the Legendre polynomials and 9 is the diffusion angle; this leads to the 
following expression for the cross-section: 

a = ^ljr(2l + l)\ ai \ 2 . (3.5) 

From the optical theorem (see Chapter 7 of [1]) we also know that the cross-section is 
proportional to the imaginary part of the amplitude in the forward direction: 

1 i fi 

a = -Im[A(0 = 0)] = — V (21 + l)|a/| 2 . (3.6) 
s s ^— ' 

1=0 

From this, we have the relation \ai\ 2 = Im(aj), that leads to the well-known condition 
(see [65]) 

|Re(a,00)| < \. (3.7) 

In the SM context, the pioneeristic work of [66] showed that, when Mjj is greater 
than a critical value ~ 1 TeV (known as unitarity bound), the elastic spherical wave 
describing the scattering of the longitudinally polarised vector bosons at very high energy 
(y/s —¥ oo) violates the condition in equation 3.7, and the perturbative stability of the 
theory breaks down. 

In the past, several efforts have been devoted to apply these methodology to a variety of 
models, in order to extract any possible information on their allowed parameter space. 
In particular, it has been already applied to scenarios with extended scalar sectors yet 
with same gauge structure as the SM, like those with additional singlets (for example, 
see [67]), doublets (for example, see [68] and [69] for non-Supersymmetric scenarios and 
[70] for Supersymmetric ones), triplets (for example, see [71]). It has also been shown 
that this approach is successful with respect to U(l) gauge group extensions of the SM 
(for example, for the case of Eq superstring- inspired minimal U(l) extensions, see [72]). 
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Then, our aim is to show how this methods could also be successfully applied to the 
minimal B — L case, taking into account the presence of two Higgs fields and four massive 
vector bosons. 

In fact, by mean of the so-called Equivalence Theorem we are allowed to compute the 
amplitude of any process with external longitudinal vector bosons Vl (V = W^, Z, Z'), 
in the high energy limit My <C s, by substituting each one of them with the related 
Goldstone bosons v = w^, z, z' and its general validity is proven (see [40]); schematically, 
if we consider a process with four longitudinal vector bosons: A(VlVl —> VlVl) = 
A(vv — > vv) + 0(m v /s). 

The intermediate vector boson exchange does not play a fundamental role in the Higgs 
boson(s) limits 6 , hence, as we intimated in Subsection 2.3.1, we simplify our approach 
by employing a theory of interacting would-be Goldstone bosons v = w , z, z' described 
by the scalar Lagrangian in Appendix A. 

In order to study the unitarity constraints in the B — L model, we calculate the tree-level 
amplitudes for all two-to-two processes involving the full set of possible (pseudo) scalar 
fields (the most relevant subset is given by table 3.1). 

Given a tree- level scattering amplitude between two spin-0 particles A(s,6), where 6 
is the scattering (polar) angle, we know that the partial wave amplitude with angular 
momentum I is given by 



It turns out that only I = (corresponding to the spherical partial wave contribution) 
leads to some bound, so we will not discuss the higher partial waves any further. 

It is well known (and we have verified) that, in the high energy limit, only the four-point 
vertices (related to the four-point functions of the interacting potential, equations (A.l)- 
(A.3) of Appendix A) contribute to the J = partial wave amplitudes, and this is 
consistent with many other aforementioned works that exploit the same methodology. 

Hence, we present the main results of our study focusing only on the relevant subset 
of all spherical partial wave amplitudes that is shown in table 3.1. Here, we should 
notice that, as one can conclude from direct computation, in the high energy limit the 
contributions in table 3.1 ticked with ~ are just a double counting of the channels ticked 
with yj or combinations of them. 

6 This is not generally true in gauge group extensions, nevertheless for this particular purpose we 
assume that g[ is perturbative and small, in such a way that any i-channel represents a subleading 
contribution to the scattering amplitude. 
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Table 3.1: The most relevant subset of two-to-two scattering processes in the minimal 
B — L model in the Higgs and would-be Goldstone boson sectors. The rows (columns) 
refer to the initial(fmal) state (or vice versa). The symbol ~ refers to processes that 
can be computed by appropriate rearrangements of those symbolised by y/. 

Moreover, the main contributions come from the so-called scattering eigenchannels, i.e., 
the diagonal elements of the "matrix" in table 3.1. In particular, for our choice of 
method, only zz — > zz and z'z' — > z'z' , and to a somewhat lesser extent also h\h\ — > h\h\ 
and h 2 h 2 — > h 2 h 2 , play a relevant role. For completeness, we list here all the arj's, 
eigenchannel by eigenchannel 7 : 

ao(zz —> zz) = 



ao(w + w — > w + w ) 



ao(z z —> z z ) 



ao(hihx — > h\hi) = 



+ 



7 Actually, in the high energy limit, ao(w + w~ — > w + w~) differs from equation (3.10) by a quantity 
~ aw due to photon and Z-boson exchange in the t-channel, but since we are applying the condition in 
equation (3.7) and aw <IC h, this correction does not change the picture of our Higgs boson mass limit 
search. 



^r- [< + K + (K - K) cos ( 2 «)] > ( 3 - 9 ) 



32M W 



aw 
16M& 



[M 2 hl + Ml + (M 2 hi - Ml) cos (2a)] , (3.10) 



327TX 2 

3aiy 



[Ml + Ml - {Ml - M 2 h2 ) cos (2a)] , (3.11) 



32M2, 



[Ml + Ml + (Ml - Ml) cos (2a)] cos 4 a 



(Ml - Ml) sin 3 (2a) 



MM W y/nx 

^ [Ml - (Ml - Ml) cos 2 a] sin 4 a, (3.12) 
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a (hih 2 ->■ hifi2) = 

+ 
+ 

a (/i2^2 — > ^2^2) = 



+ 

where 



We remark upon the fact that in the high energy limit, \/s — > oo, only the ao partial wave 
amplitude (i.e., the four-point function as one can conclude by direct comparison between 
equations (3.9)-(3.14) and equations (A.1)-(A.3) in Appendix A) does not vanish, instead 
it approaches a value depending only on M^, M^ 2 and a. Therefore, by applying the 
condition in equation (3.7), we can obtain several different (correlated) constraints on 
the Higgs masses and mixing angle, i.e., we can find the M^-M^-a subspace in which 
the perturbative unitarity of the theory is valid up to any energy scale. 

The most relevant scattering channels for the unitarity analysis are pure-z and pure 
z'-bosons scatterings. As one can see from equations (3.9)-(3.11), the limit coming from 
these two channels is unaffected by the transformation a —> —a, hence it is not restrictive 
to consider the half domain a £ [0, ^] only. Furthermore, we remind the reader that we 
are still not allowing the inversion of the Higgs mass eigenvalues, i.e., we still require 
M hx <M h2 . 

Afterwards, we analyse the space of the parameters a, M/ ll and Mh 2 , once it has been 
specified by the unitarity condition applied to the spherical partial wave scattering am- 
plitudes listed in the previous Section in the very high energy limit. 

We want to start our analysis in the M^-M^ subspace, hence we "slice" the 3-dimensional 
parameter space we are dealing with by keeping the Higgs mixing angle fixed. 



K, - (Ml - Ml) cos 2 a] sin 2 (2a) 



64-7TX 2 

[M 2 hl - {Ml - Ml) sin 2 a] sin 2 (2a), (3.13) 



3-i/ayK 



Ml - Mi) cos 2 a] cos 4 a 



UM W yfKx 
3a\y 



[Ml - mi 

M 2 hi -Ml) sin 3 (2a) 



16M 2 , 



[Ml - {Ml - Ml) sin 2 a] sin 4 a, (3.14) 



aw 



Mw 
irv 2 



(3.15) 
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Figure 3.2: Higgs bosons mass limits in the B — L model coming from the unitarity 
condition |Rc(ao)| < \ applied to the zz — > zz and z'z' — > z'z' scatterings for several 
values of x, for a = 0.01 (3.2(a)), a = 0.1 (3.2(b)), a = tt/4 (3.2(c)) and a = 0.9 tt/2 
(3.2(d)). The (blue) horizontal shadowed region corresponds to the unphysical con- 
figuration Mh ± > Mh 2 - The (red) vertical shadowed region is excluded by the LEP 

experiments. 



By applying the unitarity constraint to the spherical partial waves listed in the previous 
Section, one discovers that for a mixing angle a such that 



arctan 



M w \ 



< a < 



(3.16) 
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the allowed parameter space is completely defined by the zz — > zz eigenchannel. 

We will call "high-mixing domain" the parameter space defined by a choice of the mixing 
angle in this range, while the "low-mixing domain" is the complementary one. For 
example, since x > 3.5 TeV following the LEP analyses (Subsection 3.1.3), if we choose 
x to be exactly 3.5 TeV, then we say that the high-mixing domain, in this particular 
case, is the one for 0.073 < a < § (and, conversely, the low-mixing one is in the interval 
< a < 0.073). 

We can appreciate how the size of the Higgs mixing affects the limits on the Higgs masses 
by looking at figure 3.2, in which we plot the allowed space for the latter, limitedly to 
the two eigenchannels zz — > zz and z'z' — > z'z' , for four different values of a and three 
of x (the latter affects only the limit coming from the z'z' — > z'z' scattering). 

We see that in both expected, the light Higgs mass upper bound does not exceed 

the SM one (which is ~ 700 GeV, according to [65]), and it runs to the experimental 
lower limit from LEP (figure 3.1) as the heavy Higgs mass increases. This is because the 
two Higgses 'cooperate' in the unitarisation of the eigenchannels so that, if one Higgs 
mass tends to grow, the other one must become lighter and lighter in order to keep the 
scattering matrix elements unitarised. 

While we are in the high-mixing domain, as in figure 3.2(b)-3.2(c)-3.2(d) (where a = 0.1, 
a = j, a = 0.9 ^, respectively 8 ), the allowed region coming from the zz — > zz scattering 
is completely included in the z' z' — > z'z' allowed area, and the highest value allowed for 
the heavy Higgs mass only depends on the mixing angle via 



When we move to figure 3.2(a) (where a = 0.01, low-mixing domain) we are able to 
appreciate some interplay between the two scattering processes. In fact, in this case, 
while the zz — > zz scattering eigenchannel allows the existence of a heavy Higgs of more 
than 10 TeV, the z'z' — > z'z' scattering channel strongly limits the allowed mass region, 
with a "cut-off" on the heavy Higgs mass almost insensible to the light Higgs mass (and 
the value of the mixing angle, since we are in the low-mixing domain), that is 



For the last of these values of the mixing angle, the lower limit from LEP experiments on the light 
Higgs boson is My ll > 40 GeV, while for the first it is almost equal to the SM lower limit (Mh 1 > 115 
GeV), as illustrated in figure 3.2. 




(3.17) 




(3.18) 
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which is in agreement (under different theoretical assumptions, though) with the result 
in [72]; from a graphical point of view, in figure 3.2(a) the (green) hollow area represents 
the allowed configuration space at x = 3.5 TeV, while at x = 10 TeV the allowed portion 
of the subspace increases until the (green) double-lines shadowed region, finally 

the constraint relaxes to the (green) single line shadowed region when x = 35 TeV. 

This interplay effect arising (somewhat unintuitively) for Higgs low-mixing is due to 
the fact that the consequent decoupling between the two Higgs states requires the 
light(heavy) Higgs state to independently keep the scattering matrix elements of the 
z ( ) z { ) ^O^C) process unitary, thus realising two separate constraints: the first on 
the light (S'M-like) Higgs mass due to the zz — > zz unitarisation and the second on the 
heavy Higgs mass due to the z'z' — > z'z' unitarisation. 

To summarise, given a value of the singlet Higgs VEV x (compatible with experiment), 
the upper bound on the light Higgs boson mass varies between the SM limit and the 
experimental lower limit from LEP as long as the upper bound for the heavy Higgs 
mass increases. Moreover, when a assumes a value included in the high-mixing domain, 
the strongest bound comes from the unitarisation of the z-boson scattering, whilst in 
the low-mixing domain the bound on the heavy Higgs mass coming from that channel 
relaxes and the unitarisation induced by the z'-boson scattering becomes so important 
to also impose a cut-off (which depends linearly on x) on the heavy Higgs mass. 

This is a very important result, because it allows us to conclude that, whichever the 
Higgs mixing angle, both Higgs boson masses of the B — L model are bounded from 
above. As examples of typical values for the heavy Higgs mass, in table 3.2, we show 
some upper bounds that universally apply (i.e., no matter what the mixing angle is) 
once the singlet Higgs VEV is given. 



x (TeV) 


Max(M^) (TeV) 


3.5 


~ 10 


7 


~ 20 


10 


~ 30 


20 


~ 60 


35 


~ 100 



Table 3.2: Universal upper bound on the heavy Higgs mass, Mh 2 , in the B — L model 
as a function of the singlet Higgs VEV, x. 

Before we move on, it is also worth re-emphasising that, if the Higgs mixing angle is such 
that we are in the high-mixing case, the upper bound on the heavy Higgs boson mass 
coming from z-boson scattering is more stringent than the one coming from z'-boson 
scattering and it is totally independent from the chosen singlet Higgs VEV . 
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(c) (d) 

Figure 3.3: Heavy Higgs boson mass limits plotted against the mixing angle a in 
the minimal B — L model. We have applied the unitarity condition |Re(ao)| < h on 
zz — > zz (red straight line), z'z' — > z'z' (green dashed line) and /12/12 —> (blue 
dashed-dotted line) scatterings. This has been plotted for two hxed values of the light 
Higgs boson mass (M hl = 150 GeV (3.3(a), 3.3(c)) and M hl = 700 GeV (3.3(b), 
3.3(d))) and of the singlet Higgs VEV (x = Mz'/{2g[) = 3.5 TcV (3.3(a), 3.3(b)) and 
x = M Z '/{2g' 1 ) = 35 TeV (3.3(c), 3.3(d))). 

Nowadays, it is important to refer in our analysis to the possibility of a Higgs boson 
discovery at the Large Hadron Collider (LHC). Thus, if we suppose that a light or heavy 
Higgs mass has been already measured by an experiment it is interesting to study 
the a-Mh 2 parameter space, to see whether an hitherto unassigned Higgs state can be 
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consistent with a minimal B — L scenario. 

To this end, in figure 3.3 we fix Mh 1 and x at two extreme configurations: we take 
M/u = 150 GeV as minimum value (conservatively, taking a figure that is allowed by 
the experimental lower bound established by LEP for a SM Higgs boson) and Mh x = 700 
GeV as maximum value (close to the maximum allowed by unitarity constraints, as we 
saw before). Then, we take x = 3.5 TeV as minimum value (that is, the lower limit 
established by LEP data for the existence of a Z' of B — L origin) and M/ ll = 35 TeV as 
maximum value (that is, one order of magnitude bigger than the smallest VEV allowed 
by experiment). 

Even in this case we can separate the 2-dimensional subspace in a low-mixing region and 
a high- mixing region, as before. We can identify the first (second) as the one in which 
the upper bound is established by unitarisation through the z'(z)-boson scattering. The 
value of the mixing angle that separates the two regions in this case is given by 



Once the light Higgs boson mass is fixed, we can see how the heavy Higgs boson mass is 
bounded from above through the value defined by equation (3.18) through the z'z 1 — > z'z' 
channel, and this occurs in the low-mixing region. In particular, we can notice how the z'- 
constraining function reaches a plateau and overlaps with the ft,2^2 — > ^2^2 eigenchannel 
bound. Moreover, if we pay attention to the high-mixing region, we see that, if M/ ll is 
fixed to some low value, then the bound on the heavy Higgs mass relaxes much more as 
the mixing gets smaller and smaller with respect to the the situation in which M/ ll is 
large, where the unitarisation is shared almost equally by M^ 2 and M^. 

3.2.2 Triviality and vacuum stability bounds 

It is very well known that because of quantum corrections, the parameters which appear 
in the B — L Lagrangian (as well as in the SM Lagrangian) are energy-scale dependent 
(see [22] for details). 

This is also true for the quartic Higgs couplings as well as the gauge couplings which will 
be monotonically increasing with the energy scale, and this growth leads to constraints 
on the parameter space. 

It is important to highlight the fact that we have previously defined the minimal B — L 
by setting g = in equation (2.6). However, g is scale-dependent, and it increases with 




(3.19) 



Chapter 3. The Higgs sector parameter space 



35 



the energy-scale, so it is more appropriate to say that we define the minimal B — L by 
setting g = at the EW energy-scale. 

From this, it is clear that a correct estimate of the evolution of the couplings must 
include the g contribution. 

Therefore, the variation of the parameters is described by the Renormalisation Group 
Equation (RGE), and from [22, 73] we know that the running of Aj's is described by: 



d(Ai) 



1 



(24A 2 + A 2 + 12A 1% 2 - 9A l5 2 - 3A l5 2 - 3Ai3 2 ) , 



d(logA) 16tt 2 
d^iAO * ^(l0A| + A 2 + 4A 2 Tr[( y -) 2 ]-24A 2 (,' 1 ) 2 ) 

#4\ " ^(6A 1+ 4A 2 + 2A 3 ), 
«(log A) 8tt z 

and the running of the gauge couplings is described by: 



d(gi) l 

d(logA) 16vr 2 

rf(ffi) 1 

d(logA) 16vr 2 

d{g) 1 

d(logA) 16vr 2 



41 
~6 



.91 



12 5 f + 2-1^5 + ^9i9 
3 b 

(g 2 + 2 5 2 ) + 2^ 5 ' 1 (5 2 + 9 2 i) + Wi9 
b 3 



16 



.12- 



41 / 



~2 



(3.20) 



(3.21) 



Once we know the RGEs in equations (3.20)-(3.21), we must fix the boundary conditions 
for the evolution of the parameters. We choose the EW gauge coupling to be set by 
experiment (gi(EW) ~ 0,36) as well as the "top"-Yukawa (yt — 1). 

Regarding the neutrinos, for simplicity we consider them degenerate and we fix their 
masses to M*j~ = 200 GeV (hence, y M ~ m„ h /xV2) 9 , since this value has been proven 
to be related to interesting phenomenological scenarios (see [16]). 

Moreover, another boundary condition is set by the definition of the model: g{EW) = 0. 

Focusing on the scalar sector, the free parameters in our study are then M/^, M^ 2 , a 
and x. The general philosophy is to fix in turn some of the free parameters and scan 
over the other ones, individuating the allowed regions fulfilling the aforementioned set 
of conditions (as we previously did in Subsection 3.2.1). 



9 For semplicity we assume that yt and y M do not evolve, nevertheless their evolution could in principle 
affect quantitatively (not qualitatively) our conclusions (see [22] for further details). 
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We firstly define a parameter to be "perturbative" for values less than unity. This is a 
conservative definition, as we could relax it by an order of magnitude and still get values 
of the parameters for which the perturbative series will converge 10 . 

Then, RGE evolution can constrain the parameter space of the scalar sector in two 
complementary ways: from one side, the couplings must be perturbative, 

< Ai, 2)3 (Q0 < 1 VQ' <Q, (3.22) 

and it is usually referred to as the "triviality" condition; on the other side, the vacuum 
of the theory must be well-defined at any scale, that is, to guarantee the validity of 
eqs. (2.12) and (2.13) at any scale Q' < Q: 

< Ai j2 (Q0 and 4X 1 (Q')X 2 (Q') - A|(Q') > V Q' < Q , (3.23) 

and it is usually referred to as the "vacuum stability" condition. 

Given the simplicity of the scalar sector in the SM, the triviality and vacuum stability 
conditions can be studied independently and they both constrain the Higgs boson masses, 
providing an upper bound and a lower bound, respectively. In more complicated models 
than the one considered here, it might be more convenient to study the overall effect 
of equations (3.22)-(3.23), since there are regions of the parameter space in which the 
constraints are evaded simultaneously. This is the strategy we decided to follow. 

Figure 3.4 shows the allowed region in the parameter space M/ ll -M/ l2 for increasing 
values of the mixing angle a, for fixed VEV x = 7.5 TeV and heavy neutrino masses 
M Ufi = 200 GeV, corresponding to Yukawa couplings whose effect on the RGE running 
can be considered negligible. 

For a = 0, the allowed values for M/i 1 are the SM ones and the extended scalar sector 
is completely decoupled. The allowed space is therefore the simple direct product of 
the two, as we can see in figure 3.4(a). When there is no mixing, the bounds that we 
get for the new heavy scalar are quite loose, allowing a several TeV range for Mh 2 , still 
depending on the scale of validity of the theory. We observe no significant lower bounds 
(i.e., Mh 2 > 0.5 GeV), as the right-handed Majorana neutrino Yukawa couplings are 
negligible. 

As we increase the value for the angle, the allowed space deforms towards smaller values 
of . If for very small scales Q of validity of the theory such masses have already been 
excluded by LEP, for big enough values of Q, at a small angle as a = 0.1, the presence 

10 Notice that the parameters upon which the perturbative expansion is performed are usually of the 
form yfa — g/^/A^r, rather then being g itself. 
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Figure 3.4: Allowed values in the M hl vs. M h2 space in the B — L model by eqs. (3.22) 
and (3.23), for (3.4(a)) a = 0, (3.4(b)) a = 0.1, (3.4(c)) a = tt/4 and (3.4(d)) a = tt/3. 
Colours refer to different values of Q/GeV: blue (10 3 ), red (10 7 ), green (10 10 ), purple 
(10 15 ) and cyan (10 19 ). The shaded black region is forbidden by our convention Mh 2 > 
Mh 1 , while the shaded red region refers to the values of of the scalar masses forbidden 
by LEP. Here: x = 7.5 TcV, M Vh = 200 GcV. 

of a heavier boson allows the model to survive up to higher scales for smaller h\ masses 
if compared to the SM (in which just h\ would exist). Correspondingly, the constraints 
on Mh 2 become tighter. Moving to bigger values of the angle, the mixing between hi 
and h,2 grows up to its maximum, at a = vr/4, where h\ and hi both contain an equal 
amount of doublet and singlet scalars. The situation is therefore perfectly symmetric, as 
one can see from figure 3.4(c). Finally, in figure 3.4(d), we see that the bounds on M^ 2 
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Figure 3.5: Allowed values in the M^ 2 vs. a space in the B — L model by eqs. (3.22) 
and (3.23), for (3.5(a)) M hl = 100 GeV, (3.5(a)) M hl = 120 GeV, (3.5(a)) M hl = 160 
GeV and (3.5(a)) M hl = 180 GeV. Colours refer to different values of Q/GeV: blue 
(10 3 ), red (10 7 ), green (10 10 ), purple (10 15 ) and cyan (10 19 ). The plots already encode 
our convention Mh 2 > Mh 1 and the shaded red region refers to the values of a forbidden 
" by LEP. Here: x = 3.5 TcV, M„ h = 200 GeV. 

are getting tighter, approaching the SM ones, and those for M/ ll are relaxing. That is, 
for values of the angle 7r/4 < a < vr/2, the situation is qualitatively not changed, but 
now /12 is the SM-\ike Higgs boson. Visually, one can get the allowed regions at a given 
angle it/ 2 — a by simply taking the specular figure about the = M^ 2 line of the plot 
for the given angle a. 

Complementary to the previous study, we can now fix the light Higgs mass at specific, 
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experimentally interesting 11 , values, i.e., M/ ll = 100, 120, 160 and 180 GeV, and show 
the allowed region in the M^ 2 vs. a plane. This is done in figure 3.5. 

From this figures it is clear that the transition of hi from the new extra scalar to the 
SM-like Higgs boson as we scan on the angle. As we increase M/ ll (up to Af/ M = 160 
GeV), a bigger region in Mh 2 is allowed for the model to be valid up to the Plank scale 
(the most inner regions, in cyan). Nonetheless, such a region exists also for a value of the 
light Higgs mass excluded by LEP for the SM, M/ ll = 100 GeV, but only for big values 
of the mixing angle. No new regions (with respect to the SM) in which the model can 
survive up to the Plank scale open for > 160 GeV, as the allowed space deforms 
towards smaller values of Mh 1 . 

Combining the results of this Subsection with those on the previous one (Subsection 3.2.1), 
we are now in a position to investigate the production and decay phenomenology of both 
Higgs states of the minimal B — L model at present and future accelerators. This will 
be realised in Chapter 4. 

3.2.3 Unitarity and triviality constraints on g[ 

As well as constraining the Higgs sector parameters, the perturbative unitarity and 
-RGE'-based techniques are also helpful to constraining the g[ parameter. 

For this reason, in this Subsection, we apply this methods to the analysis of the B — L 
gauge coupling. 

Focusing on the computational aspects of this analysis, in Subsection 3.2.1 we have 
already seen that perturbative unitarity violation at high energy occurs only in vector 
and Higgs bosons elastic scatterings, our interest is focused on the corresponding sectors. 

We have also already made use of the Equivalence Theorem, which guarantees that we 
can replace the gauge bosons and Higgses interacting Lagrangian with the corresponding 
Goldstone and Higgs bosons effective theory (further details of this formalism can be 
found in Subsection 3.2.1 and in [21].) 

Moreover, since we want to focus on g' x limits, we assume that the two Higgs bosons 
of the model have masses such that no significant contribution to the spherical partial 
wave amplitude will come from the scalar four-point and three-point functions. Taking 
Higgs boson masses smaller than the unitarity limits is therefore a way to exclude any 
other source of unitarity violation different from the possible largeness of the g[ gauge 
coupling. 

n The chosen values maximise the probability for the decays hi — > bb, hi — > 77, hi — > W + W~ and 
hi — > ZZ, respectively. 
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In the search for the g[ upper limits, we assume that we can neglect all the gauge 
couplings in the covariant derivative but g[, and the equation (2.7) becomes: 

D M ~S M + 2i5i2^. (3.24) 

The inclusion of g[ in the covariant derivative gives rise to two new Feynman rules with 
respect to the set described in Appendix A, i.e.: 

Z'hiz = -2ig[ sin aip^ - p%,), (3.25) 
Z'h 2 z' = 2*0' 1 cosa(p£ | -i£), (3.26) 

where all the momenta are considered incoming and z' is the would-be-Goldstone boson 
associated with the new Z' gauge field. 

Now that the background is set, we focus on the techniques that we have used to ob- 
tain the unitarity bounds (Subsection 3.2.1) in combination with the RGE analysis 
(Subsection 3.2.2). 

Firstly, we recall that equations (3.21) and the boundary conditions where gi(EW) ~ 
0.36 and g(EW) = fully fix the evolution of g[ with the scale. 

In the search for the maximum g'^EW) allowed by theoretical constraints, we previously 
showed that the contour condition 

Si (A) < k, (3.27) 

also known as the triviality condition, is the assumption that enables to solve the system 
of equations and gives the traditional upper bound on g[ at the EW scale (where it is 
usually assumed either k = 1 or k = v4~7r, calling for a coupling that preserves the 
perturbative convergence of the theory). 

Nevertheless, what is important to say is that this is an "ad hoc" assumption, while our 
aim is to extract the boundary condition by perturbative unitarity arguments, showing 
that under certain conditions it represents a stronger constraint on the domain of g[. 

While in the search for the Higgs boson mass bound it is widely accepted to assume 
small values for the gauge couplings and large Higgs boson masses, for our purpose we 
reverse such argument with the same logic: we assume that the Higgs boson masses are 
compatible with the unitarity limits and we study the two-to-two scattering amplitudes 
of the whole scalar sector, pushing the largeness of g[ to the perturbative limit. 
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By direct computation, also in this case it turns out that only J = (corresponding 
to the spherical partial wave contribution) leads to some bound, and the only divergent 
contribution to the spherical amplitude is due to the size of the coupling g[ in the 
intermediate Z' vector boson exchange contributions. Hence, the only relevant channels 
are: z'z' — > h\hi, z'z' — > h\h2, z'z' — > h 2 h 2 - 

As an example, we evaluate the ao partial wave amplitude for z'z' — > h\h\ scattering in 
the s > M z ,,M hl limit. 

Firstly, we know that: 



4s 



M(s, cos 9) ~ (2g' 1 sin a) 2 (l - — 

\ s (1 — cos V) + IMt, cos a 



(3.28) 



and by mean of the integration proposed in equation (3.8), we extract the J = partial 
wave: 



ao(z'z' — > h\hi) 



167T 



l + 21og 



sin 2 a. 



(3.29) 



It is important to notice that the mass of the Z' acts as a natural regularisator that pre- 
serves both the amplitude and the spherical partial wave from any collinear divergence. 

Considering the three aforementioned scattering channels, their spherical partial wave 
(in the high energy limit s Mz>, Mf ll 2 ) is represented by the following matrix: 



ao = f{g'i,s;x) 



V2 







sin a 



sm a cos a 



\ 



1 2 

2 cos a 



1 • 2 

n sin a 









1 sin a cos a i cos 2 a 



where, according to equation (2.32) 
f(g[,s;x) 



(Ml! 

16vr 



l + 21og 



V2 









(2g[x) 



(3.30) 



(3.31) 



and the elements of the matrix are related to the system consisting of —^z z , -^h\hi, 
hih 2 , -j^h 2 h 2 . 

The most stringent unitarity bound on the g[ coupling is derived from the requirement 
that the magnitude of the largest eigenvalue combined with the function f(g[,s;x) does 
not exceed 1/2 (see equation (3.7)). 
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If we diagonalise the matrix in equation (3.30), we find that the maximum eigenvalue 
and the corresponding eigenvector are: 



Combining the informations of equations (3.31)-(3.32), together with the perturbative 
unitarity condition in equation (3.7), we obtain: 



In the last equation, s represents the scale of energy squared at which the scattering is 
consistent with perturbative unitarity, i.e. s = A 2 , where A is the evolution energy scale 
cut-off. 

Finally, if we consider the contour of this inequality, we find exactly the boundary 
condition that solve the set of differential equations in (3.21), giving us the upper limit 
for g[ at the EW scale. 

Evaluating the set of differential equations (3.21) by means of the well-known Runge- 
Kutta algorithm, and imposing both the unitarity and triviality condition as a two-point 
boundary value with a simple shooting method 13 , we have made a comparison between 
the two arguments for several values of x, the Higgs singlet VEV: the results are plotted 
in figure 3.6. 

By direct comparison of the two formulae in equations (3.27)-(3.33), it is easy to see 
that the unitarity bounds become more important than the triviality bounds when 



with the assumption that Mz> ~ x and k = 1. 

From this equation, it is straightforward to see that the "ad hoc" choice of the triviality 
parameter k is crucial for establishing which limit is the most stringent one. 

If we then choose a value of the VEV x that is compatible with the experimental limits 
and still in the TeV range, x £ [3.5,35] TeV (see equation (3.3)), we find that the 
unitarity bounds are more stringent than the triviality ones when the energy scale is 

12 It is important to mention the fact that the maximum eigenvalue does not depend on a, hence also 
the unitarity condition is a-independent. 

13 It consisted in varying the initial conditions in dichotomous-converging steps until the unitarity 
bound was fulfilled. 




(3.32) 




(3.33) 




(3.34) 
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Figure 3.6: Triviality (assuming k = 1, dottcd-dashcd line) and unitarity (continuous 
and dashed lines) limits on the g[ coupling of the minimal B~L model, plotted against 
the energy cut-off A in log 10 -scale, for several choices of the singlet VEV (x = 3.5 TeV, 
black lines; x = 10 TeV, red (dark grey) lines; x = 35 TeV, green (light grey) lines). 

greater than a critical value of A c ~ 10 6 GeV, and this is consistent with the results in 
figure 3.6. 

In order to summarise these results, in table (3.3) we present a comparison between 
the triviality and the unitarity bounds for several energy scales and (B — L)-breaking 
VEVs. 

Though these results are scale dependent, we see that, if A 3> A c , our method basically 
refines the triviality bound by an absolute value of ~ 0.1, that represents a correction 
of (at least) 20% on the results that have recently appeared in the literature (see [22]). 
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Log w (A/GeV) 


7 


10 


15 


19 


TB, g[(k = l) 


0.595 


0.501 


0.407 


0.357 


UB, g[(x = 3.5 TeV) 


0.487 


0.360 


0.269 


0.230 


UB, g[(x = 10 TeV) 


0.510 


0.368 


0.273 


0.232 


UB, g[(x = 35 TeV) 


0.542 


0.379 


0.277 


0.234 



Table 3.3: Triviality bounds, equation (3.27) with k = 1, and unitarity bounds, 
equation (3.33) with x = 3.5, 10, 35 TeV, for g[ in the minimal B — L model for several 

values of the energy scale A. 

In conclusion, we have shown that, by combining perturbative unitarity and RGE meth- 
ods, one can significantly constrain the g[ coupling of the minimal B — L extension of 
the SM , by imposing limits on its upper value that are more stringent than standard 
triviality bounds 14 . 



3.2.4 The fine-tuning constraint 

In this Subsection we discuss the "naturalness" aspects of the model through the theo- 
retical indication on the Higgs mass value that comes from the "fine-tuning" problem: 
this is introduced by the analysis of the one-loop radiative corrections to the Higgs boson 
mass(es). 

In the SM framework, the dominant contributions to the quantum corrections are 

3A 2 

M 2 H = (M° ) 2 + (Ml + 2M$y + M\ - 4M 2 ), (3.35) 

where A is the energy-scale cut-off and Ml is the bare mass contained in the unrenor- 
malised Lagrangian. This induced correction is quadratically divergent rather than the 
usual logarithmic ones. If A is large (for example of the order of the grand unification 
scale), a fine adjustment between the bare mass and the corrections is needed in order 
to have a physical Higgs boson with a mass in the range of EW symmetry breaking 
scale: this is the "naturalness" problem of the SM. 

However, following the Veltman conjecture (see [74]), if the following relation holds: 

Ml ~ 4M t 2 - 2M^ - Mf, (3.36) 

not only the fine-tuning problem is smoothed, but we could have an indirect indication 
of the SM Higgs mass (~ 320 GeV). 

14 Moreover, as unitarity is more constraining than triviality, the stability of the perturbative solution 
obtained through the former is already guaranteed by the latter. 



Chapter 3. The Higgs sector parameter space 



45 



We apply this argument to the B — L model: firstly, we have two Higgs bosons and two 
related expressions for the mass radiative corrections: 

A< = 

= 64Mf|U 2 (3M ^ (3M| ' + 8G?1) V) + M ' 2 (3M| ' + 8(9 ' l)2v2) 

12 (4M f 2 Mf, + w4(g[) 2 M 2 h v 2 - M%, {2M& + M 2 Z + Av 2 ((g[) 2 - 2c 2 S 2 (^) 2 ))) 

+ 4 (Ml (3M 2 , - 8(g[) 2 v 2 ) + 3 (l6c 4 a (g[) 2 M 2 h v 2 

+ Ml {-AM 2 + 2M 2 , + Mf - 4v 2 ((g[) 2 + 2c 2 > 2 (^) 2 )))) cos (2a) 

- 12( 3 ;)M Z ^ (m£ - M 2 2 + leV^M^^)) sin (2a) 

+ (M£ - M 2 2 ) ((3Mf, + 8(^)V) cos (4a) + 2(^)M^i;sin(4a))) , (3.37) 



and 



AM 2 2 = 
A 2 



[M 2 hi (3M 2 , + 8(g' 1 ) 2 v 2 ) + 3 (M 2 2 (3M|, + 8(g[ 



64M 2 ,vr 2 7; 2 

+ 4 (-16c 4 (^) 2 M 2 ^ 2 + M 2 , (-4M 2 + 2M 2 W + M 2 + 4, 2 ((^) 2 - 2c 2 S 2 (^) 2 )))) 

- 4 (M 2 2 (3M 2 , - S^) V) + 3 {l6ci(g[) 2 M 2 h v 2 

+ M 2 , (-4M 2 + 2M 2 , + M 2 - 4^ 2 ((g[) 2 + 2 C 2 S 2 (y-) 2 )))) cos (2a) 

+ l2(g[)M z ,v (-Ml + M 2 2 + l6V2c s a M uh s a v(y u )) sin (2a) 

(M 2 t - M 2 2 ) ((3M|, + S^) V) cos (4a) + 2(g[)M zl v sin (4a))) , (3.38) 

where we have chosen the neutrinos being degenerate (M^ 2 ' 3 ~ M v h, 2/1,2,3 ~ V") anc ^ 
their mixing angles being averaged (s a i, a 2, a 3 ~ s a ). 

By mean of the Veltman conjecture, we set the radiative one-loop corrections to be small 
(~ 0), obtaining a set of two equations in two variables, M^ x and M^, and the solution 
is: 



Ml — 

~ (3 {(g[)M z ,v [\2c\M 2 vh - 8M 2 + AM 2 , + 2M 2 - 3M 2 , - 16c 2 * 2 ?; V) 2 ) 
+ (gi)M Z 'V {-12ciM 2 h - 8M 2 + AM 2 , + 2M 2 + 3M 2 , - 16c 2 * 2 ?; V) 2 ) cos (2a) 
+ 8c 4 a (g[) 2 M 2 h v 2 sin (2a) + M 2 , (-4M 2 + 2M 2 V + M 2 - 2(g[) 2 v 2 ) sin (2a) 
8c 2 M|, S y(^) 2 sin(2a) - A^clM vh s a v(y u ) (2 cos (a) sec (2a) sin (a) x 
x (8(31) V cos 2 (a) + 3Mf, sin 2 (a)) + ^{g'^Mz'V sin ( 2a ) tan (2a)))) -r 
-j- (-12(j;)M Z( d cos (2a) + (-3M|, + 4^) V) sin (2a)) , (3.39) 
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and 




~ (3 {(g'^Mz'V {~l2c A a Ml h + 8M 2 - 4M 2 , - 2M 2 + 3Mf, + 16c 2 S 2 ?; V) 2 ) 

+ (j^M^i; {-Ylc\Ml h - 8M 2 + 4M 2 , + 2M 2 + 3M|, - 16c 2 a s 2 a v 2 (y») 2 ) cos (2a) 

+ 8 c 4 a (g[) 2 M 2 h v 2 sin (2a) + Aff, {-AM 2 + 2M 2 , + M§ - 2(0!) V) sin (2a) 

- &c 2 a Ml,s 2 a v 2 (y u ) 2 sin (2a) + A^f2c\M vh s a v(y v ) (2 cos (a) sec (2a) sin (a) x 

x (3M|,cos 2 (a) + 8(c/ 1 )Vsin 2 (a)) - 3(ffi)M Z r?j sin (2a) tan (2a)))) -r- 

-7- (— 12^) Myu cos (2a) + (-3Mf, + 4(si) V) sin (2a)) . (3.40) 



Indeed, it is important to say that the Veltman conjecture is only indicative of the one- 
loop behaviour 15 , and evaluating the stability of the Veltman solution at two-loops is 
not among the purposes of this work. However, it is interesting to notice that in the 
decoupled scenario a = 0, the solution is: 



We find the SM-like solution plus a non-significative contribution from heavy neutrinos 
(equation (3.41)), while the second solution is totally ruled by the interplay between the 
heavy neutrinos and Z' boson mass (equation (3.42)). 

Again, too many free-parameters are floating around, hence we can only give an illus- 
trative example of a choice of parameters related to a reasonable amount of fine-tuning. 

The role of heavy neutrinos in the first equation is absolutely negligible, so as in the SM 
case M/ ll ~ 320 GeV represents a good choice for the light Higgs; assuming a relatively 
heavy right-handed neutrinos (~ 400 GeV) and a relatively light Z' boson (~ 700 GeV) , 
the Veltman conjecture implies that a theory with a heavy Higgs boson with mass of 
~ 470 GeV is the most "natural" . 




~ 4M 2 -2M^-M 2 +8 C 2 ^V) 2 



(3.41) 




a 9 3ME, 



(3.42) 



15 Moreover, since too many free parameters enter in the game, no significative study could be per- 
formed in the B — L context. 



Chapter 4 

Higgs phenomenology at colliders 



This Chapter is devoted to the presentation of the results of our phenomenological 
investigation of the B — L Higgs sector. 

Firstly, after a summary of our choice of parameters and computational details, we 
present the total decay width and the branching ratios (BRs) of the Higgs bosons, in 
order to set the background. 

Then, we focus on the LHC experimental scenario, with emphasis on the possible peculiar 
signatures of the model. 

Finally, we will present a set of possible signature cross sections that could represent a 
relevant test for the SM break-down at future LCs. 



4.1 Parameter space and implementation of the model in 
CalcHEP 

As spelled out in the previous Chapter, the independent physical parameters of the 
Higgs sector of the scenario considered here are: 

• M/i i; Mh 2 and a, the Higgs boson masses and mixing angle. We investigate this 
parameter space spanning over continuous intervals in the case of the first two 
quantities while adopting discrete values for the third one (taking into account all 
the experimental and theoretical constraints presented in Chapter 3). 

In order to explore efficiently the expanse of parameter space pertaining to the minimal 
B — L model, we introduce two extreme conditions, which makes the model intuitive, 
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though at the end it should be borne in mind that intermediate solutions are phenomeno- 
logically favoured. The two conditions are obtained by setting: 

1. a = 0, this is the "decoupling limit", with h\ behaving like the SM Higgs. 

2. a = 7£, this is the so-called "inversion limit", with h% behaving like the SM 
Higgs (though recall that this possibility is phenomenologically not viable, see 
Subsection 3.1.1). 

Furthermore, concerning the strength of Higgs interactions, some of the salient phe- 
nomenological behaviours can be summarised as follows 1 : 

• SM-like interactions scale with cosa(sina) for h]_(h,2)', 

• those involving the other new B — L fields, like Z' and heavy neutrinos, scale with 
the complementary angle, i.e., with sina(cosa) for h\{h2)\ 

• triple (and quadruple) Higgs couplings are possible and can induce resonant be- 
haviours, so that, e.g., the hi — > h\ h\ decay can become dominant if M^ 2 ~ 2Mf ll . 

Other than Mh 1 , M^ 2 and a, additional parameters are the following: 

• g[, the new U{\)b-l gauge coupling. We will adopt discrete perturbative values 
for this quantity, taking into account the analysis of the allowed parameter space 
of Subsection 3.2.3. 

• Mz', the new gauge boson mass, considering the experimental limits presented in 
Subsection 3.1.3 (equation (3.3)). 

• M Vh , the heavy neutrino masses, are assumed to be degenerate, diagonal and 
relatively light. 

• M Vv the SM (or light) neutrino masses; they have been conservatively taken to 
be M Vl = 10~ 2 eV in order to fulfil the cosmological bound requirement presented 
in Subsection 3.1.3; they are assumed to be degenerate and diagonal. 

The numerical analysis was performed with CalcHEP [37] with the model introduced 
through LanHEP [34]. Moreover, the implementation has been enriched with the fol- 
lowing plug-ins: 

All of these results are a consequence of the way the particles couple, i.e., the Feynman rules in 
Appendix B. 
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• the one-loop vertices g — g — /ii(/i2) 2 , 7 — 7 — ^1(^2) and 7 — Z(Z') — h±(h2) via W 
gauge bosons and heavy quarks (top, bottom and charm) have been implemented, 
adapting the formulae in [77]. 

• Running masses for top, bottom and charm quarks, evaluated at the Higgs bo- 
son mass: Q = M^ 1 (M/ l2 ), depending on which scalar boson is involved in the 
interaction. 

• Running of the QCD coupling constant, at two-loops with 5 active flavours. 

4.2 Branching ratios and total widths 

Moving to the Higgs boson decays, figure 4.1 shows the BRs for both the Higgs bosons, 
h\ and h2, respectively. Only the two-body decay channels are shown here. 

Regarding the light Higgs boson, the only new particle it can decay into is the heavy 
neutrino (we consider a very light Z' boson unlikely), if the channel is kinematically 
open. In figure 4.1(a) we show this case, for a small heavy neutrino mass, i.e., M Uh = 50 
GeV, and we see that the relative BR of this channel can be rather important, as 
the decay into 6-quark pairs or into W boson pairs, in the range of masses 110 GeV 
< M/ ll < 150 GeV. Such range happens to be critical in the SM since here the SM 
Higgs boson passes from decaying dominantly into 6-quark pairs to a region in masses 
in which the decay into W boson pairs is the prevailing one. These two decay channels 
have completely different signatures and discovery methods. The fact that the signal 
of the Higgs boson decaying into 6-quark pairs is many orders of magnitude below the 
natural QCD background spoils its sensitivity. In the case of the B — L model, the 
decay into heavy neutrino pairs is therefore phenomenologically very important, besides 
being an interesting feature of the B — L model if M Uh < Myy , as it allows multileptons 
signatures of the light Higgs boson. Among them, there is the decay of the Higgs boson 
into 31, 2j andfJ T (that we have already studied for the Z' case in reference [16]), into 
4£ &ndfJ T (as, again, already studied for the Z' case in reference [18]) or into M and 2j 
(as already studied, when I = (j,, in the 4 th family extension of the SM [78]). All these 
peculiar signatures allow the Higgs boson signal to be studied in channels much cleaner 
than the decay into 6-quark pairs. 

In the case of the heavy Higgs boson, further decay channels are possible in the B — L 
model, if kinematically open. The heavy Higgs boson can decay in pairs of the light 
Higgs boson (hy ~~ ^ h±h\) or even in triplets (/12 — > hih\hi), in pairs of heavy neutrinos 

2 Finally, the NLO QCD fc-factor for the gluon fusion process [2, 75, 76] has been used. Regarding the 
other processes, we decided to not implement their fc-factors since they are much smaller in comparison. 
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FIGURE 4.1: (4.1(a)) Branching ratios for hi for a = 2tt/5 and M Uh = 50 GeV 
and (4.1(c)) hi total width for a choice of mixing angles and (4.1(b)) BRs for 
for a = 3tt/20 and M hl = 120 GeV, M z > = 210 GeV and M Vh = 200 GeV and (4.1(d)) 
hi total width for a choice of mixing angles. 



and Z' bosons. Even for a small value of the angle, figure 4.1(b) shows that the decay of 
a heavy Higgs boson into pairs of the light one can be quite sizeable, at the level of the 
decay into SM Z bosons for = 120 GeV. It is important to note that this channel 
does not have a simple dependence on the mixing angle a, as we can see in figure 4.2. 

The BRs of the heavy Higgs boson decaying into Z' boson pairs and heavy neutrino 
pairs decrease as the mixing angle increases, getting to their maxima (comparable to 
the W and Z ones) for a vanishing a, for which the production cross section is however 
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FIGURE 4.2: Dependence on the mixing angle a of (4.2(a)) BR(h 2 — > h\ hi), of (4.2(b)) 
BR(h 2 -> v h v h ) and of (4.2(c)) BR(h 2 -> Z' Z'). 



negligible. As usual, and also clear from figure 4.1(b), the decay of the heavy Higgs 
boson into Z' gauge bosons is always bigger than the decay into pairs of fermions (the 
heavy neutrinos, even when summed over the generations as plotted), when they have 
comparable masses (here, Mz> = 210 GeV and M Uh = 200 GeV). 

The other standard decays of both the light and the heavy Higgs bosons are not modified 
substantially in the B — L model (i.e., the decay of the Higgs boson to W boson pair 
is always dominant when kinematically open, while before that the decay into 6-quarks 
is the prevailing one; further, radiative decays, such as Higgs boson decays into pairs of 
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Figure 4.3: Dependence on the mixing angle a of the three body decays (4.3(a)) 
BR(h 2 -> h x hi hi) and (4.3(b)) BR(h 2 -> h x V V) (V = W ± ,Z) for M hl = 120 GeV 
and (4.3(c)) for M hl = 90 GeV, respectively. 



photons, peak at around 120 GeV). Only when other new channels open, the standard 
decay channels alter accordingly. This rather common picture could be altered when 
the mixing angle a approaches vr/2, but such situation is phenomenologically not viable 
(see Subsection 3.1.1). 

Figures 4.1(c) and 4.1(d) show the total widths for h\ and h%, respectively. In the first 
case, few thresholds are clearly recognisable, as the heavy neutrino one at 100 GeV 
(for angles very close to 7r/2 only), the WW and the ZZ ones. Over the mass range 
considered (90 GeV < M hl < 250 GeV), the particle's width is very small until the WW 
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threshold, less than 1 — 10 MeV, rising steeply to few GeV for higher h\ masses and 
small angles (i.e., for a S'M-like light Higgs boson). As we increase the mixing angle, 
the couplings of the light Higgs boson to SM particles is reduced, and so its total width. 

On the contrary, as we increase a, the hi total width increases, as clear from figure 4.1(d). 
Also in this case, few thresholds are recognisable, as the usual WW and ZZ gauge boson 
ones, the light Higgs boson one (at 240 GeV) and the ti one (only for big angles, i.e., 
when hi is the S'M-like Higgs boson). When the mixing angle is small, the hi total 
width stays below 1 GeV all the way up to M^ 2 ~ 300 -f- 500 GeV, rising as the mass 
increases towards values for which Th 2 ~ Mh 2 ~ 1 TeV and hi loses the meaning of 
resonant state, only for angles very close to n/2. Instead, if the angle is small, i.e., less 
than 7r/10, the ratio of width over mass is less than 10% and the heavy Higgs boson is 
a well defined particle. In the decoupling regime, i.e., when a = 0, the only particles hi 
couples to are the Z' and the heavy neutrinos. The width is therefore dominated by the 
decay into them and is tiny, as clear from figure 4.1(d). 

As already mentioned, figure 4.2 shows the dependence on the mixing angle a of the BRs 
of hi into pairs of non-SM particles. In particular, we consider the decays hi — > hi h\ 
(for two different h\ masses, = 90 GeV and Mj ll = 120 GeV, only for the allowed 
values of a), hi — > Vh v h an d hi — > Z' Z' (not influenced by M/ ll ). As discussed in 
Section 4.1, the interaction of the heavy Higgs boson with SM (or non-SM) particles has 
an overall sin a (or cos a, respectively) dependence. Nonetheless, the BRs in figure 4.2 
depend also on the total width, that for a > 7r/4 is dominated by the h<i — > W + W~ 
decay. Hence, when the angle assumes big values, the angle dependence of the hi 
BRs into heavy neutrino pairs and into Z' boson pairs follow a simple cot a behaviour. 
Regarding hi — > hi hi, its BR is complicated by the fact that the contribution of this 
process to the total width is not negligible when the mixing angle is small, i.e., a < 7r/4. 
In general, this channel vanishes when a — > 0, and it gets to its maximum, of around 
10% -T- 30% of the total width, as a takes a non-trivial value, being almost constant with 
the angle if it is small enough. 

The heavy Higgs boson can be relatively massive and the tree-level three-body decays 
are interesting decay modes too. Besides being clear BSM signatures, they are crucial to 
test the theory beyond the observation of any scalar particle: its self-interactions and the 
quartic interactions with the vector bosons could be tested directly in these decay modes. 
In the B — L model with no Z — Z' mixing, the quartic interactions that can be tested as 
hi decay modes, if the respective channels are kinematically open, are: hi — > hi hi hi, 
hi — > hi W + W~ and hi —> hi Z Z, as shown in figure 4.3, again for = 90 GeV and 
120 GeV. Although possible, hi — > hi Z' Z' is negligible always, even if the Z' boson is 
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light enough to allow the decay. For M z > = 210 GeV, BR(h 2 -»• h x Z' Z') < 1(T 5 for 
M h2 < 2 TeV. 

The BRs for both the h 2 — > h\ h\ h\ and the h 2 — > hi V V (V = W^, Z) channels are 
maximised roughly when the mixing between the two scalars is maximum, i.e., when a ~ 
7r/4, regardless of Mh x . The former channel, that is interesting because would produce 
three light Higgs bosons simultaneously, can contribute at most at 10 -3 of the total width 
for h 2 , as we are neglecting values of M^ 2 and a for which Yh 2 ~ M^ 2 (see figure 4.1(d)). 
For instance, for M/ l2 = 800 GeV, a needs to be less the n/5 to have a reasonable small 
width-over-mass ratio (~ 10%), and BR(h 2 — > h\h\hi) < 0.6 • 10~ 3 . The situation is 
similar for the latter channel, involving pairs of SM gauge bosons. Again, for Mh 2 = 800 
GeV and a = vr/5, BR(h 2 hi W + W~) = 2BR{h 2 -*h x ZZ) = 10~ 3 for m hl = 120 
GeV. For = 90 GeV, the mixing angle is constrained to be bigger than 7-7r/20. For 
these values and the same M^ 2 as before, such BRs are doubled. 

4.3 Higgs bosons at the LHC 

In this Section we present our results for the analysis of the scalar sector of the minimal 
B — L model at the LHC. We shortly introduce the scheduled working plan at the 
accelerator. Then, we present cross sections at -^/s = 7 and 14 TeV for the two Higgs 
bosons. Finally, we will focus on some phenomenologically viable signatures and their 
event rates. 

4.3.1 The LHC scheduled working plan 

The scheduled programme is planned to be the following: 

• 7 TeV is total energy of the two proton beams (energy in the hadronic center of 
mass) and 1 fb _1 is the scheduled integrated luminosity (~ 1 — 2 operational years 
in the time scale). This is what we label as "early discovery scenario". 

• 14 TeV is total energy of the two proton beams (energy in the hadronic center of 
mass) and 300 fb _1 is the scheduled integrated luminosity (~ 10 operational years 
in the time scale). This is what we label as "full luminosity scenario". 

4.3.2 Standard production mechanisms 

In figure 4.4 we present the cross sections for the most relevant production mechanisms, 
i.e., the usual SM processes such as gluon-gluon fusion, vector-boson fusion, ti associated 
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FIGURE 4.4: Cross sections in the minimal B — L model for h\ at the LHC (4.4(a)) 
at V« = 7 TeV and (4.4(c)) at ^ = 14 TeV, and for h 2 (4.4(b)) at ^ = 7 TcV and 
(4.4(d)) at y/s = 14 TcV. Dashed lines in figures 4.4(a) and 4.4(c) refer to a = 0. The 
dotted part of the lines in figure 4.4(d) refer to h 2 masses excluded by Unitarity (see 

Subsection 3.2.1). 

production and Higgs-strahlung. For reference, we show in dashed lines the SM case 
(only for hi), that corresponds to a = 0. 

Comparing figure 4.4(c) to figure 4.4(a), there is a factor two enhancement passing from 
a y/s = 7 TeV to a yfs = 14 TeV centre-of-mass energy at the LHC. 



The cross sections are a smooth function of the mixing angle a, so, as expected, every 
sub-channel has a cross section that scales with cos a (sin a), respectively for h\ (fo)- 
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As a general rule, the cross section for h\ at an angle a is equal to that one of I12 for 
7r/2 — a. In particular, the maximum cross section for /12 (i.e., when a = 7r/2) coincides 
with the cross section of h\ for a = 0. 

We notice that these results are in agreement with the ones that have been discussed 
in [79-81] in the context of a scalar singlet extension of the SM, having the latter the 
same Higgs production phenomenology. Moreover, as already showed in [79], also in 
the minimal B — L context an high value of the mixing angle could lead to important 
consequences for Higgs boson discovery at the LHC: h\ production could be suppressed 
below an observable rate at y/s = 7 TeV and heavy Higgs boson production could be 
favoured, with peculiar final states clearly beyond the SM, or even hide the production 
of both (if no more than 1 fb _1 of data is accumulated). Instead, at y/s = 14 TeV we 
expect that at least one Higgs boson will be observed, either the light one or the heavy 
one, or indeed both, thus shedding light on the scalar sector of the B — L extension of 
the SM discussed in this work. 

4.3.3 Non-standard production mechanisms 

All the new particles in the B — L model interact with the scalar sector, so novel 
production mechanisms can arise considering the exchange of new intermediate particles. 
Among the new production mechanisms, the associated production of the scalar boson 
with the Z' boson and the decay of a heavy neutrino into a Higgs boson are certainly 
the most promising, depending on the specific masses. Notice also that the viable 
parameter space, that allows a Higgs mass lighter than the SM limit of 114.4 GeV for 
certain a — M^ 2 configurations, enables us to investigate also production mechanisms 
that in the SM are subleading, as the associated production of a Higgs boson with a 
photon. Figures 4.5 and 4.6 show the cross sections for the non-standard production 
mechanisms, for y/s = 14 TeV and several values of a. 

Figures 4.5(a) and 4.5(b) show the cross sections for associated production with the Z' 
boson of h\ and of hi, respectively, for several combinations of Z' boson masses and g[ 
couplings. The process is 



and it is dominated by the Z' boson's production cross sections (see [16, 20]). Although 
never dominant (always below 1 fb), this channel is the only viable mechanism to produce 
/12 in the decoupling scenario, i.e., a = 0. 




(4.1) 



In figure 4.6 we plot the cross sections of the other non-standard production mechanisms 
against the light Higgs mass, for several choices of parameters (as explicitly indicated in 
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Figure 4.5: Cross sections in the minimal B — L model for the associated production 
with the Z' B _ L boson (4.5(a)) of hi at a = tt/A and (4.5(b)) of fi2 at a = 0. 

the labels). We superimposed the red-shadowed region in order to avoid any value of the 
cross section that has been already excluded by LEP constraints (see Subsection 3.1.2), 
mapping each value of the boundary cross section as produced by the related maximum 
value allowed for the light Higgs mass M/ ll (at fixed mixing angle a). 

First of the showed plots is the decay of a heavy neutrino into a Higgs boson. The whole 
process chain is 

qq^ Z' ->• u h v h ->• VhVihxp) , (4.2) 

and it requires to pair produce heavy neutrinos, again via the Z' boson (see [16, 82] for 
a detailed analysis of the pp — > Z' — > v^v^ process and other aspects of Z' and heavy 
neutrinos phenomenology in the minimal B — L model). Although rather involved, this 
mechanism has the advantage that the whole decay chain can be of on-shell particles, 
besides the peculiar final state of a Higgs boson and a heavy neutrino. For a choice 
of the parameters that roughly maximises this mechanism (Mz> = 900 GeV, g[ = 0.13 
and M Vh = 200 GeV), figure 4.6(a) shows that the cross sections for the production of 
the light Higgs boson (when only one generation of heavy neutrinos is considered) are 
above 10 fb for < 130 GeV (and small values of a), dropping steeply when the light 
Higgs boson mass approaches the kinematical limit for the heavy neutrino to decay into 
it. Assuming the transformation a —> n/2 — a, the production of the heavy Higgs boson 
via this mechanism shows analogous features. 
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Figure 4.6: Cross sections in the minimal B — L model for the associated production 
of h\ (4.6(a)) with one heavy and one light neutrinos, (4.6(b)) with a photon via 7, Z 
and Z' bosons exchange (same legend as in figure 4.6(a) applies here) and (4.6(c)) in 
the vector-boson fusion, all at s/s = 14 TeV. The red shading is the region excluded by 

LEP constraints (Subsection 3.1.2). 

Next, figures 4.6(b) and 4.6(c) shows the associated production of the light Higgs boson 
with a photon. The processes are, respectively, 



qq^j/Z/Z' -^ 1 h 1 



(4.3) 



Chapter 4. Higgs phenomenology at colliders 



59 



via the SM neutral gauge bosons (7 and Z) and the new Z' boson, and 

qq 7 /ti q" q" , (4.4) 
through vector-boson fusion (only W and Z bosons). 

In the first instance, we notice that the Z' sub-channel in equation (4.3) is always 
negligible, as there is no Z' — W — W interaction and the V — h — 7 effective vertex is 
only via a t-quark loop (an order of magnitude lower than the V — h — 7 effective vertex 
via a W boson loop) [77]. What is relevant in these two channels is that the light Higgs 
boson mass can be considerably smaller than the LEP limit (they are valid for the SM, 
or equivalently when a = in the B — L model). Hence, the phase space factor can 
enhance the mechanism of equation (4.3) for small masses, up to the level of 1 fb for 
M^j < 60 GeV (and suitable values of the mixing angle a, depending on the experimental 
and theoretical limits, see [21, 22] for a complete tratement of the allowed parameter 
space of the Higgs sector of the minimal B — L model). Moreover, it has recently 
been observed that the associated production with a photon in the vector-boson fusion 
channel could be useful for low Higgs boson masses in order to trigger events in which 
the Higgs boson decays into 6-quark pairs [83]. Complementary to that, the process in 
equation (4.3) can also be of similar interest, with the advantage that the photon will 
always be back-to-back relative to the 6-quark pair. For comparison, figures 4.6(b) and 
4.6(c) show the cross section for these processes 3 . Certainly, for a h± boson heavier than 
the SM limit, vector-boson fusion is the dominant process for associated production of 
hi with a photon, and this is also true for > 60 GeV. However, for light Higgs 
boson masses lower than 60 GeV, the two mechanisms of eqs. (4.3) and (4.4) become 
equally competitive, up to the level of 0(1) fb each, for suitable values of the mixing 
angle a. 

4.3.4 Event Rates 

In this Section we combine the results from the Higgs boson cross sections and those 
from the BR analysis in order to perform a detailed study of typical event rates for some 
Higgs signatures which are specific to the B — L model. 

Before all else, we remind the terminology previously introduced in Subsection 4.3.1: 
we will generally refer to an "early discovery scenario" by considering an energy in the 
hadronic CM of yfs = 7 TeV and an integrated luminosity of f L = 1 fb -1 (according 
to the official schedule, this is what is expected to be collected after the first couple of 

3 In order to produce figure 4.6(c), we included the following cuts: P?'^ > 15 GeV, |?7 7 | < 3 and 
|7^ e ^| < 5.5, where "jet" refers to the actual final state, though we use partons here to emulate it [83]. 
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Figure 4.7: Cross section times BR contour plot for the B—L process pp —> hi —t Vh^h 
at the LHC with y/s = 7 TeV, plotted against M hl -a, with M„ h = 50 GeV (4.7(a)) and 
M Vh = 60 GeV (4.7(b)). Several values of cross section times BR have been considered: 
(7 = 5 fb (black line), a = 10 fb (red line), a = 50 fb (green line), a = 100 fb (blue 
line) and a = 250 fb (violet line). The red-shadowed region is excluded by the LEP 

experiments. 



years of LHC running) and to a "full luminosity scenario" by considering an energy in 
the hadronic CM of \fs = 14 TeV and an integrated luminosity of f L = 300 fb -1 
(according to the official schedule, this is what is expected to be realistically collected 
at the higher energy stage). 

As we shall see by combining the production cross sections and the decay BRs already 
presented, the two different scenarios open different possibilities for the detection of 
peculiar signatures of the model: in the "early discovery scenario" there is a clear possi- 
bility to detect a light Higgs state yielding heavy neutrino pairs while the "full luminosity 
scenario" affords the possibility of numerous discovery mechanisms (in addition to the 
previous mechanism, for the heavy Higgs state one also has decays of the latter into Z' 
boson and light Higgs boson pairs). 

Firstly, we focus on the "early discovery scenario": in this experimental configuration, 
the most important B — L distinctive process is represented by heavy neutrino pair 
production via a light Higgs boson, through the channel pp — > hi — > Vh v h- In figure 4.7 
we show the explicit results for the pp — > h\ — > Vh u h process at the LHC with y/s = 7 
TeV, for M Vh = 50 GeV (figure 4.7(a)) and M Uh = 60 GeV (figure 4.7(b)), obtained 
by combining the light Higgs boson production cross section via gluon-gluon fusion only 
(since it represents the main contribution) and the BR of the light Higgs boson to heavy 
neutrino pairs. The obtained rate is projected in the M/ ll -a plane and several values 



Chapter 4. Higgs phenomenology at colliders 



61 




^ n/2 



m M = 120 GeV 

Sg— a = 400 fb 
||— a = 250 fb ■ 
It- a = 100 lb 
a = 25fb 




200 300 



500 700 1000 2000 
m h2 (GeV) 

(a) 



^1-6 

£ 

-1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 





^ n/2 



%/3 



31/4 



» jc/6 



200 300 



500 700 1000 2000 
m h2 (GeV) 

00 



FIGURE 4.8: Cross section times BR contour plot for the B—L process pp —> hi — > h\h\ 
at the LHC with yfs = 14 TeV, plotted against M h2 -a, with M hl = 120 GeV (4.8(a)) 
and Mf ll = 240 GeV (4.8(b)). Several values of cross section times BR have been 
considered: a = 10 fb (violet line), a = 25 fb (light-blue line), a = 100 fb (blue line), 
a = 250 fb (green line) and a = 400 fb (red line) . The red-shadowed region is excluded 

by unitarity constraints. 



of the cross section times BR have been considered: a = 5, 10, 50, 100 and 250 fb. 
The red-shadowed region takes into account the exclusion limits established by the LEP 
experiments. 

Even considering a low-luminosity scenario (i.e., / L ~ 1 fb" 1 ), there is a noticeable 
allowed parameter space for which the rate of such events is considerably large: in 
the case of M Uh = 50 GeV, when the integrated luminosity reaches f L = 1 fb -1 , we 
estimated a collection of ~ 10 heavy neutrino pairs from the light Higgs boson production 
and decay for 100 GeV< < 170 GeV and 0.057T < a < 0.487T, that scales up to 
~ 10 2 events for 110 GeV< M hl < 155 GeV and 0.16vr < a < 0.46vr. In the case 
of M Ufi = 60 GeV, we estimated a collection of ~ 10 heavy neutrino pairs from Higgs 
production for 120 GeV< M hl < 170 GeV and 0.067T < a < 0.487T, that scales up to 
~ 10 2 events for 125 GeV< M hl < 150 GeV and 0.25vr < a < 0.44vr. 

If we consider instead the "full luminosity scenario" , there are several important distinc- 
tive signatures: pp — > hi — > h\h\, pp — > hi — > Z'Z' and pp — > hi — > Vh u h- I R figure 4.8 
we show the results for light Higgs boson pair production from heavy Higgs boson decays 
at the LHC with y/s = 14, TeV for M hl = 120 GeV (figure 4.8(a)) and M hl = 240 GeV 
(figure 4.8(b)). Again, if we project the rates on the bi-dimensional M^-a plane, we 
can select the contours that relate the cross section times BR to some peculiar values. 
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Figure 4.9: Cross section times BR contour plot for the B—L process pp I12 —> Z'Z' 
at the LHC with ^/i = 14 TcV, plotted against M h2 -a, with M z > = 210 GeV (4.9(a)) 
and Mz< = 280 GeV (4.9(b)). Several values of cross section times BR have been 
considered: a = 0.085 fb (light-blue line), a = 0.15 fb (blue line), a = 0.3 fb (green line), 
a = 0.85 fb (red line). The red-shadowed region is excluded by unitarity constraints. 



Considering an integrated luminosity of 300 fb -1 , we can relate a = 25(250) fb to 
7500(75000) events, hence for both choices of the light Higgs mass the a-M^ 2 parameter 
space offers an abundant portion in which the event rate is noticeable for light Higgs 
boson pair production from heavy Higgs boson decays: when Mh 1 = 120 GeV the 
process is accessible almost over the entire parameter space, with a cross section peak 
of 400 fb in the 240 GeV< M h2 < 400 GeV and 0.13vr < a < 0.30vr intervals, while in 
the M/ ll = 240 GeV case the significant parameter space is still large, even if slightly 
decreased, with a cross section peak of 25 fb in the 480 GeV< Mh 2 < 800 GeV and 
0.067T < a < 0.327T region. 

In figure 4.9 we show the results for Z' boson pair production from heavy Higgs boson 
decays at the LHC with y/s = U TeV for M z > = 210 GeV (figure 4.9(a)) and M z > = 280 
GeV (figure 4.9(b)). Again, if we project the rates on the bi-dimensional M^ 2 -a plane, 
we can select the contours that relate the cross section times BR to some peculiar values. 
Here, we have that a = 0.085(0.85) fb corresponds to 25(250) events, hence for both 
choices of Z' mass the a-M^ 2 parameter space offers an abundant portion in which the 
event rate could be interesting for Z' boson pair production from heavy Higgs boson 
decays: for M z > = 210 GeV the process has a peak of 0.85 fb in the 420 GeV< M h2 < 650 
GeV and 0.037T < a < 0.257T region, while if M z > = 280 GeV a noticeable parameter 
space is still potentially accessible with a rate peak of 0.3 fb (100 events) in the 560 
GeV< M h2 < 800 GeV and 0.03vr < a < 0.19vr region. 
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Figure 4.10: Cross section times BR contour plot for the B — L process pp — >• h,2 ~ > 
Vh^h a t the LHC with yfs — 14 TeV, plotted against Mh 2 -a, with M Vh = 150 GeV 
(4.10(a)) and M Vh = 200 GeV (4.10(b)). Several values of cross section times BR have 
been considered: a = 0.085 fb (light-blue line), a = 0.15 fb (blue line), a = 0.3 fb 
(green line), a = 0.85 fb (red line). The red-shadowed region is excluded by unitarity 

constraints. 



In analogy with the previous two cases, in figure 4.10 we show the results for heavy neu- 
trino pair production at the LHC with y/s = 14 TeV plus M Uh = 150 GeV (figure 4.10(a)) 
and M Uh = 200 GeV (figure 4.10(b)). The usual contour plot displays a sizable event 
rate in the a-Mh 2 parameter space for both choices of the Vh mass. For example, when 
M Uh = 150 GeV we find a cross section times BR peak of 0.85 fb (~ 250 events) in the 
320 GeV< M h2 < 520 GeV and 0.03vr < a < 0.33tt region, while if M Ufi = 200 GeV we 
find a peak of 0.85 fb in the 450 GeV< M h2 < 550 GeV and 0.03vr < a < 0.21tt region. 



4.4 Higgs bosons at future linear colliders 



In this Section we present our results for the scalar sector of the minimal B — L model 
at future LCs. In general, sub- TeV CM energies (y/s = 500 GeV) will be suitable for 
an ILC, multi-TeV CM energies (y/s = 3 TeV) will be appropriate for CLIC while the 
case 1 TeV may be appropriate to both. In all cases, results will be shown for some 
discrete choices of the Z' mass and of the scalar mixing angle a. Their values have been 
chosen in each plot to highlight some relevant phenomenological aspects. Concerning 
single Higgs production, we will distinguish the standard production mechanisms (via 
SM gauge bosons) from the novel mechanisms present in the model under discussion 
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(emphasising in particular the role of the Z' gauge boson). Finally, plots for double 
Higgs production will also be presented. 

4.4.1 The future linear collider running proposals 

Although there are not official approvals of either ILC or CLIC yet, we already know 
what are the energy parameters of the two proposed machines: 

• ILC: 500 GeV is the planned initial energy, with the open possibility of an upgrade 
to 1 TeV, and 500 fb _1 of integrated luminosity at fixed energy; moreover, the 
possibility is planned to span over the energy range up the maximum energy (500 
GeV or 1 TeV) at 10 fb _1 of integrated luminosity. 

• CLIC: 3 TeV is the planned initial energy, with the open possibility of an upgrade 
to 5 TeV; there is no scheduled integrated luminosity yet, then we conservatively 
assume that it will correspond to the ILC prototype, i.e. 500 fb _1 . 

The standard Initial State Radiation (ISR) functions are implemented, according to the 
formulae in [19, 36, 84]. 

Finally, though the beam-strahlung parameters have only been set for the ILC prototype, 
in this work we assume that the same set of values holds for the CLIC framework, hence 
we will take into account these values as they appear in the "ILC Reference Design 
Report" (see [85]). We list them in table 4.1. 





Nominal value 


Unit 


Bunch population 


2 


xl0 lu 


RMS bunch lenght 


300 


/im. 


RMS horizontal beam size 


640 


nm 


RMS vertical beam size 


5.7 


nm 



Table 4.1: Nominal values of beam parameters at the ILC (see [85]). 
4.4.2 Standard single-Higgs production mechanisms 

Figure 4.11 shows the cross sections for the standard production mechanisms of a single 
Higgs boson: vector boson fusion (V = W ,Z), the strahlung from the Z boson and 
the associated production with a t-quark pair. These standard production mechanisms 
are modulated by a cosa(sina) prefactor in the vertices when considering h\{h,2), as 
generally true for all the scalar singlet extensions of the SM. Hence, we will not spend 
too much time in discussing them. 
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FIGURE 4.11: Cross sections for the process standard Higgs boson production mecha- 
nism as a function of the mass at the LC for a = 37r/10 for (4.11(a)) hi and at y/s = 500 
GeV, for (4.11(b)) hi at y/s = 1 TeV and for (4.11(c)) h x at s/s = 3 TeV. The dashed 

lines refer to a = 0. 

As well known [2] , for low energies, Higgs-strahlung from a Z boson is competitive with 
W-boson fusion, that becomes the main mechanism as we increase the CM energy. At 
\fs = 1 TeV, the cross sections for a light Higgs boson vary from O(0. 1) pb for M/ tl = 100 
GeV to O(10) fb for = 600 GeV, in the decoupling regime (i.e., for a = 0, that 
corresponds to recovering the SM Higgs boson). All the cross sections are then scaled 
by a factor cos a (sin a) when considering hi (h^). At this value of CM energy then, the 
associated production with a t-quark pair reaches its highest value (for the CM energies 
we are plotting), i.e., O(l) fb for M/ ll < 200 GeV, depending on the mixing angle. 
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Not surprisingly, as we increase the CM energy towards y/s = 3 TeV, VF-boson fusion 
increases considerably, staying around fractions of pb for several masses and angles, for 
both Higgs bosons, while the Higgs-strahlung from the Z boson plunges towards cross 
sections of the order of few fb. 

The associated production with a i-quark pair in the SM scenario is the least effective 
production mechanism, with cross sections of few fb at most. However, we will show 
in the following Subsection that this mode can be enhanced by the presence of the Z' 
boson. 

4.4.3 Non-standard single-Higgs production mechanisms 

In this Section we will discuss the novel mechanisms to produce a Higgs boson (both 
the light one or the heavy one) in the minimal B — L model. All the new features arise 
from having a Z' that interacts with both the scalar and fermion sectors. We recall here 
another important feature: the Z 1 boson is dominantly coupled to leptons [16]. In fact: 

^BR(Z'^U) ~ I; (4.5) 

e 

Y,BR{Z'^ q q) ~ 1; (4.6) 
q 

and in particular, BR{Z' — > e + e~) ~ 15%, which makes a lepton collider the most 
suitable environment for testing this model. 

Again, the Z' mass and g[ gauge coupling values have been chosen to respect the con- 
straints coming from LEP and Tevatron (Subsection 3.1.3). 

We start by showing the cross section for the associated production of a Higgs boson 
and a Z' boson, as in figure 4.12. Due to the stringent bounds on the Z' boson mass 
and coupling to fermions, a sub- TeV CM energy collider is not capable of benefiting 
from this production mechanism, especially because of the naive kinematic limitation 
in the final state phase space. In other words, there is not enough energy to produce a 
Z' and a Higgs boson, if both are on-shell. This is clear in figures 4.12(a) and 4.12(b), 
where a light Z' boson (with mass of 500 GeV) gives cross sections below 0.1 fb. For 
a Z' boson of 700 GeV mass instead, the cross sections can be of the order of few fb, 
only for Higgs masses below 300 GeV, the kinematical limit. These results are similar 
to the LHC, in which the Z' strahlung process has cross sections below 1 fb (below 0.1 
fb for M/j 2 > 400 GeV) [24] , even if the LC is expected to accumulate roughly an order 
of magnitude more integrated luminosity than the LHC. 
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FIGURE 4.12: Cross sections for the process e + e — » Z * — » hu^Z' (4.12(a)) for hi 
and (4.12(b)) for /i 2 at the LC at y/s = 1 TeV and (4.12(c)) for hi and (4.12(d)) for h 2 

at the LC at -Js = 3 TeV. 



The situation is considerably improved for a multi-TeV collider, not anymore limited in 
kinematics. As shown in figures 4.12(c) and 4.12(d), a Higgs boson can be produced 
in association with a Z' boson of 1.5 TeV mass with cross sections of ~ 10 fb in the 
whole range of the scalar masses considered, rising to 0(100) fb if Mz> = 2.1 TeV is 
considered (and a suitable value for g[ is chosen). Although such configuration could 
suffer from kinematical limitations for the Higgs boson mass to be produced, the cross 
sections when the scalar mass is close to 700 GeV (the maximum value considered here) 
are still above those when a Z' boson of 1.5 TeV mass is considered, regardless of the 
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Higgs boson. It is crucial to note that this is the only production mechanism that can 
potentially lead to the discovery of the heavy Higgs boson in the decoupling limit, i.e., 
for a — > 0. As previously stated, the strahlung from Z' is not suitable for the LHC, 
making a multi-TeV LC possibly the ultimate chance for its discovery. 

As anticipated, the associated production with a t-quark pair can be enhanced exploit- 
ing the Z' boson. In figures 4.13(a) and 4.13(b) are shown the cross sections for the 
associated production of a Higgs boson and a pair of i-quarks, for M// = 500 GeV, 700 
GeV and in the case of a much heavier Z' boson, hence decoupled, for y/s = 1 TeV. As 
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known, the Higgs boson in this channel can be radiated both by a t-(anti)quark or by 
the vector boson, even though the fraction of events with a Z boson is negligible with 
respect to the t-quark pair produced by a photon. Therefore, in the SM, the measure of 
the Higgs coupling to the t-quark is possible, though difficult because of the small cross 
sections [86]. We are therefore left to evaluate the relative contribution of the Z' boson, 
to check whether the same situation holds. 

It is first interesting to note that, in the decoupling limit (i.e., for vanishing scalar mixing 
angle a), h\ does not couple directly to the Z' boson. Nonetheless, the Z' boson can 
decay to t-quark pair, one of which then radiates the light Higgs boson. This channel 
has the same final state than the SM ones, and will therefore increase the total number 
of events, as clear from figure 4.13(a). Hence, the chances of measuring the (SM-like) 
Higgs boson to t-quark coupling are improved in this case, only slightly for y/s = 1 TeV 
but quite considerably for y/s = 3 TeV and a few TeV Z' boson mass. 

As we increase the scalar mixing angle, the relative contribution of the Z' boson in- 
creases, although the total cross sections for y/s = 1 TeV fall below the fraction of fb, 
making it even harder to be observed. The situation is opposite for the multi-TeV CM 
energy case (figures 4.13(c) and 4.13(d)), in which the Z' boson is produced abundantly 
and it can enhance the Higgs boson associated production with a t-quark pair. In this 
case, anyway, it is not true anymore that the majority of the events are those in which 
the Higgs boson is radiated by a t-quark: the Higgs-strahlung from the Z' boson is now 
an important channel, as clear from figures 4.12(c) and 4.12(d) and from the fact that, 
for low Z 1 masses, the total cross section is smaller as we start increasing the angle 
(due to the reduced coupling to the t-quark), while for TeV Z' boson masses it always 
increases. If the Z' boson mass is below the maximum CM energy of the collider, the 
fraction of strahlung events off the Z' boson can be reduced by tuning the CM energy 
and sitting at the peak of the Z' boson itself. In this case, the vast majority of Z' bosons 
are produced on-shell, enhancing the total cross sections and the portion of events in 
which the Z' boson decays into a t-quark pair, one of which will then radiate the Higgs 
boson. The possibility of sitting at the peak of the Z' boson is therefore very important 
phenomenologically, as it allows the Higgs coupling to the t-quark to be measured much 
more precisely than in the SM, as the cross section in the minimal B — L model for 
this channel can rise up to 10 -j- 100 fb, depending on the Higgs boson mass and mixing 
angle. Notice that for hi the angle has to be small (i.e., less than 7r/5) to allow for this 
measurement, as only in this case hi couples more to the t-quark than to the Z' boson, 
though Z' strahlung events are still important (the ratio of the two subchannels is in fact 
< 5%). This situation is exactly specular when the heavy Higgs boson is considered: 
when the CM energy is maximal, the associated production with a t-quark pair has 
good cross sections but it does not allow a direct measurement of the Higgs boson to 
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Figure 4.14: Cross sections for the associated production of the light higgs boson and 
one heavy and one light first generation neutrinos (via Z' —> Vh^h) a * the LC (4.14(a)) 
for y/s = 1 TeV and (4.14(b)) for y/s = M z >. 



t-quark coupling, that instead is possible for big angles (a > 37r/10) when sitting at the 
Z' boson peak. Notice that, in this configuration, the total cross section is independent 
of the Z' boson mass, if the t-pair and the Higgs boson can be produced all on-shell. 
Otherwise, the cross sections are suppressed by the phase space. 

Next, a possibility already highlighted in [16, 24] for the LHC is using the heavy neutrino 
as a source of light Higgs bosons. Besides to provide a further production mechanism 
and being a very peculiar feature of the minimal B — L model, it also allows for a direct 
measure for the Higgs boson to heavy neutrinos coupling when the decay of the Higgs 
to neutrino pairs is kinematically forbidden. Though, in [24] is showed that it gives low 
cross sections at the LHC, making it hard to probe. In contrast, a LC is a more suitable 
environment to test this mechanism. One reason is that the Z' couples dominantly to 
leptons, as already intimated. Further, the possibility of tuning the CM energy and 
sitting exactly on the Z' peak will enhance the Z' production cross section by a factor 
of roughly 10 3 . Another key factor is that the BR of a heavy neutrino into a light Higgs 
boson and a light neutrino is ~ 20% (at the very most, see [16]), when kinematically 
allowed. This mechanism is not suitable for the heavy scalar though: since it is heavier 
than the light one, for sure one would observe the latter first. Altogether, for a Z' boson 
of 700 GeV mass, figure 4.14(a) shows the cross sections for the production of a (first 
generation only) heavy neutrino pair and the subsequent decay of one of them into a 
light Higgs boson, for two different masses of the heavy neutrino, at \fs = 1 TeV. At 
this stage, the mechanism is giving 0(1 4- 10) fb cross sections for a heavy neutrino of 
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Figure 4.15: Cross section for e + e~ -> Z(Z')* -> (£ = e,/x). (4.15(a)) Black 

line is for M z > = 420 GeV, other lines for M z > = 700, 1500, 7000 GeV; = 500 GeV. 
(4.15(b)) Black line is for M z , = 1050 GeV, other lines for M z > = 1400, 3500 GeV; 

v/i = 1 TeV. 

200 GeV mass, decreasing to 0(1) fb when a mass of 300 GeV is considered, for a good 
range in the mixing angle. Figure 4.14(b) shows the full potentiality of this model at 
a LC: by sitting on the Z' peak, the heavy neutrino pair production is enhanced by a 
factor ~ 10 3 , giving cross sections well above the pb range for a large portion of the 
allowed parameter space, and staying above 10 fb whatever the mixing angle. When 
kinematically allowed though, this peculiar mechanism really carries the hallmark of the 
minimal B — L model and it does not depend dramatically on the Z' mass, if below the 
maximum CM energy of the collider. 

Finally, the interference between the Z and the Z' could play an important role in the 
scalar sector, besides the f-quark Yukawa coupling measure presented above. As well 
known, and remarked upon in [19], the negative interference between the neutral gauge 
bosons can be substantial. One could then look for information about a further neutral 
vector boson also by looking at the interference when a Higgs boson is radiated from 
the Z bosons. To highlight this effect, in this model it is possible to select just the 
leptonic decay modes of the vector bosons, reducing the predominance of the Z boson. 
Nonetheless, as shown in figure 4.15, such effects are minimal when the Z 1 boson mass 
is above the CM energy of the LC. 

Other subleading processes for Higgs boson production consist of single production in 
association with two vector bosons and of double Higgs boson production. In these 



Chapter 4. Higgs phenomenology at colliders 



72 



contexts, many new mechanisms arise in this model, especially exploiting the Z' boson 
and because a resonant h-2 — > h\h\ process is allowed. 

4.4.4 Single scalar production in association with a pair of vector 
bosons 

The SM gauge boson pair production has really large cross sections at the LC, therefore 
the radiation of a Higgs boson could still have observable rates. Also, once the Higgs 
boson has been seen in the main production mechanisms of Section 4.4.2, these subpro- 
cesses could be useful to test the quartic coupling to the SM gauge bosons. Figure 4.16 
shows the case for y/s = 0.5, 1 and 3 TeV CM energies. We neglect here final state 
photons: although the cross section of channels comprising the photon could be com- 
parable to the WW subchannel (see, e.g., [2]), the absence of a direct coupling to the 
scalar bosons ensures that these channels would not provide further informations than 
the other processes that have been considered. 

We then see that the WW channel is roughly an order of magnitude higher than the 
ZZ one and that, for low Higgs boson masses, the cross sections decrease as we increase 
the CM energy. However, a larger CM energy allows the production of more massive 
scalars and to avoid kinematic limitations. So that, if y/s = 1 TeV is preferable to test 
these mechanisms for Higgs boson masses between 100 and 300 GeV (with comparable 
or higher cross sections to the case for y/s = 500 GeV), the y/s = 3 TeV configuration 
is essential for masses above 300 GeV, for both the light and heavy Higgs boson. Unless 
very high (and disfavoured) values of the mixing angle, the h\WW channel has cross 
section above 0.1 fb for the whole range of scalar masses considered. The opposite is 
true for the heavy Higgs boson, for which only big values of a allow this channel to be 
above 0.1 fb. The case for the ZZ channel is different: since its cross sections are rather 
small, it has chances of being detected, staying above 0.1 fb, only for = 500 GeV 
and y/s = 1 TeV. For y/s = 3 TeV its observation requires very high statistics. 

The cross sections for the case of one Z' boson in the final state could be important 
and comparable to the WW channel. Also, this particular channel is useful to test the 
absence of a tree-level h — Z — Z' coupling. Figure 4.17 shows the cross sections for 
y/s = 3 TeV for two values of the Z' mass, Mz> = 1-4 and 2.1 TeV, and suitable g[ 
coupling. The heavier the Z' boson, the higher the cross sections, until kinematical 
limitations occur. In fact, the cross sections for Mz> = 2.1 TeV are always above those 
for Mz' = 1.4 TeV for scalar masses below 600 GeV, for which the process with the 
lighter Z' boson overtakes. It is important to note that the behaviour of these processes 
with the scalar mixing angle is opposite to the previous case. Hence, for hi and for 
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Figure 4.16: Cross sections for the light Higgs boson production with two vector 
bosons (V = W , Z) as a function of the mass at the LC for a = 37r/10 (4.16(a)) for 
h! at = 500 GeV, (4.16(b)) for ht at = 1 TeV and (4.16(c)) for h x at = 3 
TeV. The dashed lines refer to a = 0. 



small values of the angle, the associated production with a pair of SM gauge bosons is 
favoured, while the process with the Z' boson is favoured for big angles. For /12 it is 
again the opposite. For intermediate angles, instead, both processes can have small but 
observable rates, between 0.1 and 0(1) fb, for both Higgs bosons and in the whole range 
of considered masses. 



Finally, notice that the case for a = represents the radiative correction to the Z(Z') 
strahlung of h\{h2) for a Z'{Z) boson emission. 
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Figure 4.17: Cross sections for the Higgs boson production with a Z and a Z' boson 
as a function of the scalar mass at the LC (4.17(a)) for hi and (4.17(b)) for /12, at 
\/s = 3 TeV. The dashed lines refer to a = 0. 



4.4.5 Double scalar production mechanisms 



The observation of a Higgs boson pair is crucial to measure parameters of the scalar 
Lagrangian entering directly in the trilinear and quartic self-couplings [87-89] , although 
it requires high statistics and large CM energy. Remarkable in this sense is the possible 
complementarity between the LHC and LCs, as shown in [90, 91]. 

In the minimal B — L model, the /12 — > h\h\ process is also present for a large portion of 
the parameter space, contrary to the MSSM case, for instance (where it is important 
only for very low values of tan/3, region that has been constrained at LEP [92]). On 
the one side, light Higgs boson pair production is enhanced by this channel, especially 
when it is resonant. On the other side, also the hi — h\ — h\ coupling is directly testable. 
Moreover, the Z' boson can give further scope to this, providing an extra mechanism 
for Higgs pair production, both without and through heavy Higgs boson production. 

Figure 4.18 shows the standard production mechanisms of a pair of light Higgs bosons. 
In the SM case (or when we neglect /12), they are the same mechanisms discussed 
in Subsection 4.4.2 when a further Higgs boson is attached. Cross sections for these 
processes are always below 0.1(1) fb at -^/s = 1(3) TeV, and above 0.1 fb at -^/s = 3 
TeV only for the H^-fusion process and for M/ ll < 350 GeV, as clear from figures 4.18(a) 
and 4.18(c), respectively. 
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Figure 4.18: Cross sections for the double light higgs boson production at the LC 
(4.18(a)) alone and (4.18(b)) via h 2 , for sfs~ = 1 TeV and (4.18(c)) alone and (4.18(d)) 
via hi, for \/s = 3 TeV. The dashed lines in figures 4.18(a) and 4.18(c) refer to a = 0, 
while in figures 4.18(b) and 4.18(d) refer to a — ir/20. 



When instead the light Higgs boson pair is originated by the decay of the heavy Higgs 
boson, cross sections can be of 0(1) — O(10) fb at y/s = 1 — 3 TeV, in the VF-fusion 
channel (at most, when the mixing is maximal, i.e., a ~ 7r/4, and we chose M^ 2 = 500 
GeV.). If we choose a higher value for M^ 2 , a more massive hi boson can be pair 
produced, but with a sensibly lower cross sections: for Mh 2 = 700 GeV, they are roughly 
a factor 5(2) smaller than for M/j 2 = 500 GeV at y/s = 1(3) TeV. Notice that the cross 
sections are constant with Mj ll as long as the hi —> h\h\ decay is allowed. This is a 
consequence of having chosen a specific value for M^ 2 and that BRQ12 — > hihi) ~ 20% 
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FIGURE 4.19: Cross sections for the process e + e~ -> H 2 Z' HiHiZ' (4.19(a)) for 
•s/s = 1 TeV and (4.19(c)) for ^/s = 3 TeV, for suitable values of Mh 2 and for the 
process e + e~ — > H\H\Z' (4.19(b)) at y 7 ? = 3 TeV, several values of the angle and of 

M z >. 



is approximately constant for M/ ll > M\y, Mz [24]. Other channels give smaller cross 
sections: the Z-fusion channel is always an order of magnitude below the VF-fusion one, 
while the strahlung from the Z boson channel gives at most ~ 0.2(0.08) fb for M^ 2 = 500 
GeV and ~ 0.08(0.06) fb for M h2 = 700 GeV at ^/s = 1(3) TeV, respectively. 

As anticipated, the Z' boson in the minimal B — L model can give further scope to 
produce also a pair of light Higgs bosons at a LC, both directly (without or through /12) 
and indirectly (pair producing heavy neutrinos). 
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Figure 4.19 shows double Higgs strahlung from the Z' boson (for y/s = 3 TeV only) and 
the case in which hi is radiated from the Z' boson and it subsequently decays into a light 
Higgs boson pair. The double Higgs-strahlung from the Z' boson at y/s = 1 TeV has 
negligible cross sections, below 10 -3 fb, especially because of kinematic limitations, and 
therefore we neglect it here. The cross sections for double Higgs-strahlung at y/s = 3 TeV 
are presented in figure 4.19(b), where we see that, for Mz> = 2.1 TeV (and g[ = 0.3), 
a pair of light Higgs bosons can be produced with cross section > 0.1 fb for < 300 
GeV and for big values of the scalar mixing angle (roughly bigger than vr/4). The 
situation improves if we consider the Higgs-strahlung of hi from the Z' boson and its 
subsequent decay into h\ pairs. Notice that this channel reduces 4 as we increase the 
value of the mixing angle, vanishing in the decoupling regimes (both for a = and 
7r/2). At y/s = 1 TeV this process is still limited by the kinematics: the higher the 
Z' boson mass the higher the cross sections and the smaller the producible hi mass. 
For Mz> = 700 GeV (and suitable values for the g[ coupling), the light Higgs boson 
can be pair produced through hi with cross sections bigger than 0.1 fb (for a < vr/4, 
up to 4 fb) through a heavy Higgs boson of 250 GeV, hence for hi masses up to 120 
GeV only. To extend the range in M/ il , a higher mass for the heavy Higgs boson has to 
be considered, needing a smaller Z' boson mass: the cross sections in this case become 
unobservable, below 10 -2 fb. If the collider CM energy is increased though, heavier 
h^s can be pair produced through the heavy Higgs boson, in association with a much 
heavier Z' boson, with bigger cross sections. Figure 4.19(c) shows that, for Mz' = 2.1 
TeV, a heavy Higgs boson with 500 GeV mass can pair produce the light Higgs boson 
with cross sections well above the fb level up to M/ ll = 200 GeV, reaching O(10) fb for 
small (but not negligible) values of the mixing angle (i.e., 7r/20 < a < vr/5). If a Z' 
boson of 1.5 TeV mass is considered, there are no more kinematical limitations for the 
producible hi boson and, in the case of M^ 2 = 700 GeV, an even heavier h\ can be pair 
produced, up to masses of 350 GeV with cross sections bigger than 0.1 fb and 0(1) fb 
for small (but not negligible) values of the mixing angle (i.e., for the same values of the 
previous case). 

The high cross sections of figure 4.14 (and the fact that BR{vh — > h±ui) ~ 20%) allows 
one to consider the case in which both heavy neutrinos decay into a light Higgs boson 
each. In figure 4.20 we show this case. Once again, the possibility of tuning the CM 
energy of the LC to sit at the Z' boson peak is crucial to test this mechanism. Without 
it, the cross sections would be about 0.1 4- 1 fb for small values of the scalar mixing 
angle only: for instance, for a = 67r/20, light Higgs boson pair production through 
heavy neutrino pair production (via the Z' boson, one generation only) are above 0.1 fb 

4 This is true when both the Z' boson and the heavy Higgs boson are on-shell. When /12 is an off-shell 
intermediate state, the cross sections for light Higgs pair production via hi increases as we increase the 
value of the mixing angle. 
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Figure 4.20: Cross sections for the associated production of two light higgs bosons 
and two light first generation neutrinos (via Z' — > Vh^h) at the LC (4.20(a)) for y/s = 1 

TeV and (4.20(b)) for y/s= M z , . 



for M/ij < 125 GeV, for both heavy neutrino masses chosen. When instead the CM is 
tuned to the Z' boson peak, the cross sections are enhanced and well above the fb level 
whatever the value for the mixing angle, for h\ masses kinematically allowed, reaching 
a few pb (or fractions of pb) for small mixing angles and values. 

Finally, the two Higgs bosons could be produced together, as shown in figure 4.21. 
Although subleading, this mechanism is peculiar for several reasons: it requires both 
Higgs bosons to be (simultaneously) significantly coupled to the gauge bosons, it has 
a very complicated dependence upon a, due to the trilinear and quartic scalar self- 
couplings, that makes it not invariant under a — > 5 — a (being a the scalar mixing 
angle) 5 and it is maximum when the mixing is maximal, i.e., for a = 7r/4. These 
processes could be important to reconstruct the scalar potential and the whole set of 
self- interaction couplings, together with the /12 — > h\h\ decay. 

If at \fs = 1 TeV the cross sections for this process are always below 0.1 fb, at y/s = 3 
TeV the Fd^-fusion channel can produce the two Higgs bosons with cross sections of 
fractions of fb, up to 0.08(0.02)^-0.2(0.3) fb, for M hl < M fl2 = 300(500) GeV. The double 
Higgs-strahlung from a Z' boson of 2.1 TeV mass (and g[ = 0.3) has also comparable 
cross sections, of O(0. 1) fb for < M^ 2 = 300 GeV only, for values of the mixing angle 
close to maximal. Notice that the cross sections for this process scale approximately 

5 The behaviour of the trilinear and quartic self-interaction couplings with the scalar mixing angle 
is not a simple trigonometric function, as one can verify by checking the Feynman rules listed in Ap- 
pendix B. 
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FIGURE 4.21: Cross sections for the associated production of the two higgs bosons at 
the LC through the standard production mechanisms (4.21(a)) for Mh 2 = 300 GeV 
and (4.21(b)) for Mh 2 = 500 GeV, for y/s = 3 TeV, and in association with a Z' boson 
(4.21(c)) for M h2 = 300 GeV and (4.21(d)) for M h2 = 500 GeV, for ^/s = 3 TeV and 

several Z' masses. 



with sin 2a, whatever production mechanism is considered. The mixing angle can be 
measured from other processes and used as an input for these channels, provided that 
also both scalar masses have been measured elsewhere. If so, the deviation of the cross 
sections from the naive ones (when the two Higgs bosons are produced independently, i.e., 
neglecting the self-interactions, that would be exactly proportional to sin 2a) will give 
further indications about the self-interaction couplings. Very high statistics is required 
for such a study, barely within the potentiality of the next generation of LCs. 
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Conclusions 



In this Thesis we have investigated the phenomenology of the Higgs sector of the minimal 
B — L model at present and future colliders. 

The model is realised by mean of a minimal extension of the SM, i.e., by gauging the 
broken 1 U{1)b-l symmetry in addiction to the SM gauge symmetry 2 , hence realising 
the following symmetry pattern: 

SU(3) C x SU(2) L x U{1) Y x 17(1) B -l- (5-1) 

This extension is triple-minimal: in the gauge boson sector (adding one gauge boson, 
Z'), in the Higgs sector (adding one scalar boson, /12) and in the fermion sector (adding 
one heavy right-handed neutrino per generation). 

In Chapter 1 we have presented the main motivations that call for the U(1)b-l exten- 
sion, in particular the phenomenological necessity of giving mass to the light neutrinos 
(via the "see-saw" mechanism) through a TeV-scale symmetry breaking. We have also 
established two possible experimental environments in which this extension could be 
studied: the LHC (a currently operative collider) and the ILC/CLIC (two future linear 
collider proposal). 

In Chapter 2 we have presented the formal aspects of the minimal B — L model, focusing 
on the details that are needed to implement the model Lagrangian in any of the public 
softwares that are devoted to Feynman Rules generation (such as LanHEP, FeynRules, 
etc.). 

1 A± the TeV energy scale. 

2 Assuming no-mixing between the two gauge groups at tree-level. 



81 



Chapter 5. Conclusions 



82 



In Chapter 3 (based on the results appeared in [21-23]) we have presented a detailed 
study of the Higgs sector parameter space. In particular, we have focused on the analysis 
of the parameter space, with emphasis on the experimental and theoretical exclusion 
methods, in order to set the basis for the phenomenological analysis of the Higgs sector 
at colliders. 

As for the theoretical methods, we have presented a full analysis on unitarity bounds in 
the Higgs sector of the minimal B — L model. Using the equivalence theorem, we have 
evaluated the spherical partial wave amplitude of all possible two-to-two scatterings in 
the scalar Lagrangian at an infinite energy, identifying the zz — > zz and z'z' — > z'z' 
processes as the most relevant scattering channels for this analysis (z^ is the would-be 
Goldstone boson of the vector boson). Then, we have shown that these two channels 
impose an upper bound on the two Higgs masses: the light one cannot exceed the SM 
bound while the limit on the heavy one is established by the singlet Higgs VEV, whose 
value is presently constrained by LEP and could shortly be extracted by experiment 
following a possible discovery of a Z' . We also studied how the discovery of a light Higgs 
boson at the LHC could impact on the heavy Higgs mass bounds in the minimal B — L 
model and we discovered that the lighter the h\ mass the more loose is the bound on 
M/j 2 , except in the low-mixing region (a — > 0) of the Higgs parameter space, in which 
the knowledge of the x VEV is again fundamental. 

Then, we have investigated the triviality and vacuum stability conditions of the minimal 
B — L model with a particular view to define the phenomenologically viable regions of 
the parameter space of the scalar sector, by computing all relevant RGEs at the one- 
loop level in presence of all available experimental constraints. The RGE dependence 
on the Higgs masses and couplings (including mixings) has been studied in detail for 
discrete choices of the singlet Higgs field VEV, in order to make a fruitful comparison 
with the unitarity case. 

Thereafter, we have shown that, by combining perturbative unitarity and RGE methods, 
one can significantly constrain the g[ coupling of the minimal B — L extension of the SM, 
by imposing limits on its upper value that are more stringent than standard triviality 
bounds. 

Finally, we have made an illustrative study of the "fine-tuning" induced by the quan- 
tum one-loop corrections of the Higgs boson masses, giving the correct solution for the 
Veltman conjecture in the minimal B — L context. 

In Chapter 4 (based on the results appeared in [24, 25]), we have studied in detail the 
phenomenology of the Higgs sector of the minimal B — L model at colliders. 
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After a short discussion about the implementation of the model in CalcHEP, we have 
presented the Higgs bosons branching ratios and total widths for some interesting points 
of the parameter space. 

Then, we have investigated both the foreseen energy stages of the LHC (and correspond- 
ing luminosities). While virtually all relevant production and decay processes of the two 
Higgs states of the model have been investigated, we have eventually paid particular at- 
tention to those that are peculiar to the described B — L scenario. The phenomenological 
analysis has been carried out in presence of all available theoretical and experimental 
constraints and by exploiting numerical programs at the parton level. While many Higgs 
signatures already existing in the SM could be replicated in the case of its B — L ver- 
sion, in either of the two Higgs states of the latter (depending on their mixing), it is 
more important to notice that several novel Higgs processes could act as hallmarks of 
the minimal B — L model. These include Higgs production via gluon-gluon fusion, in 
either the light or heavy Higgs state, the former produced at the lower energy stage of 
the CERN collider and decaying in two heavy neutrinos and the latter produced at the 
higher energy stage of such a machine and decaying not only in heavy neutrino pairs 
but also in Z' and light Higgs ones. For each of these signatures we have in fact found 
parameter space regions where the event rates are sizable and potentially amenable to 
discovery. Our results have laid the basis for the phenomenological exploitation of the 
Higgs sector of the minimal B — L model at the LHC. 

Finally, we have studied the potential of future LCs in establishing the structure of the 
Higgs sector of the minimal B — L model. We have considered both an ILC and CLIC. 
The scope of either machine in this respect is substantial as a large variety of Higgs pro- 
duction processes are accessible. The latter include both single and double Higgs boson 
channels, at times produced in association with heavy particles, both SM (W and Z 
bosons and t (anti)quarks) and B — L ones (Z' boson and Vh neutrinos), thus eventually 
yielding very peculiar signatures at detector level. This variety of accessible Higgs pro- 
duction processes potentially allows future LCs to accurately pin down the structure of 
the B — L Higgs sector, including not only the masses and couplings of both Higgs states 
pertaining to this scenario, but also trilinear and quartic self-couplings between the two 
scalar bosons themselves. On this score, the interplay and complementarity of measures 
at the LHC and at LCs is fundamental. The extension in the gauge sector, with a Z' 
boson dominantly coupled to leptons, is fundamental to distinguish this model from the 
classic scalar extensions of the SM. Although the scalar Lagrangian is rather a simple 
one, we showed that the new signatures and production mechanisms led by the Z' are 
quite peculiar and not shared with any extension of the SM that keeps its gauge content 
minimal. Finally, the fermion sector can also have very important consequences for the 
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B — L scalar sector discovery and identification, allowing for peculiar Higgs bosons decay 
patterns. 

In conclusion, we firmly believe that the analysis that we have presented in this Thesis 
could lay the basis for the phenomenological exploitation of the Higgs sector of the 
minimal B — L model at present and future colliders, representing a "today's challenge" 
for the former's working schedule and a pressing motivation for the latter's approval. 

For the sake of completeness, it is of fundamental importance to point out the fact 
that a significative number of open issues is still under investigation: in the first place, 
one of the next step of our research is to study the phenomenology of the minimal 
B — L at future photon-photon colliders. Secondly, a whole set of next-to- leading- 
order corrections should be evaluated both at the LHC and LCs. Thereafter, it should 
be important to investigate the impact of the mixing between the two U{1) groups, 
realising the non-minimal B — L model. Finally, one of the most important open issues 
lies in the supersymmetrisation of the model, for the purpose of making another step of 
our bottom-up approach, toward the inclusion of the minimal B — L model in a bigger 
theoretical grand unification picture. 



Appendix A 



The scalar potential 



In this Appendix, we rewrite the interaction part of equation (2.10) in terms of mass 
eigenstates, separating four-point and three-point functions and classifying them by the 
nature of the involved fields. 

The part of the interacting potential that contains four-point functions involving only 
would-be Goldstone bosons is: 
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The part of the interacting potential that contains four-point functions involving both 
would-be Goldstone and Higgs bosons is: 
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The part of the interacting potential that contains four-point functions involving only 
Higgs bosons is: 
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The part of the interacting potential that contains three-point functions involving both 
would-be Goldstone and Higgs bosons is: 
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The part of the interacting potential that contains three-point functions involving only 
Higgs bosons is: 
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Appendix B 

The minimal B 



— L Feynman rules 



In this Appendix we list the Feynman rules of the minimal B — L model in the Feynman- 
gauge; the labelling of the fields is straightforward, and follows the notation that has 
been introduced in Chapter 2. 

We remark upon the fact that the following list of Feynman rules has been generated 
by means of the LanHEP package. 

All the vertices must be coupled to a phase "i" . 

All the momenta appearing in the vertices are incoming. 



• e is the electric charge. 

• s w {c w ) => sin 9w (cos Qw)- 

• s a (c a ) =>• sina(cosa). 

• V is the CKM matrix (see [93, 94]). 

• s a i(c a i) is the sinus(cosinus) of the "see-saw" mixing of the i th neutrino generation 
(no mixing between generations has been considered). 
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Hi 




z v 


c a eM\Y gUV 

C w Sw 


Hi 




z 


1 ic a e ( » v tf) 


Hi 


z z 




— — ( 4s*,,Cn,Mu/ A1 — Sn,eXix) 


Hi 




Z' v 


-8s a <7i 2 x 5 ^ 


Hi 




z' 




Hi 


z' z' 


2— ( S^cXoX — S-,^C^M\aiX'x\ 

e \ ^ ^ a -t va yy s\j y 


H 2 


H 2 


H 2 


— 3^ ( 4s,„s 3 Mm/Ai + 2c 3 eXr>x + s 2 c^eA^x 
+2c 2 s„,Sn,MwA'i') 


H 2 


w+ 


w- 




H 2 


w + 

A* 






H 2 


w + 


W~ 




H 2 


w + 




-i(4s tl) s a il%Ai + c a eA 3 2;) 
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Fields in the vertex 



Variational derivative of Lagrangian by fields 



H 2 Z v 
H 2 z 

H 2 z z 

H 2 Z n Z v 

H 2 Z\ z' 

H 2 z' z' 

fa Tb Hi 

fa T b H 2 

fa Tb Z/j, 

fa T b Z 

fa T b Z'n 

fa Ul\ W~ 

f a Vl\ W~ 

fa Vh\ W~ 

f a vh\ W~ 

Pa M6 \ 

Pa fJ-b Hi 

Pa M& H 2 

Pa fJ-b Zfj, 

P-a M6 Z 

Pa fJ-b Z'n 

Pa Vl\ W~ 

pa Vl\ W~ 

Pa Vh\ W~ 



— 2 Q 



1 ies 



2 Cin Si 



(Ps - Pi) 



-\{As w s a M w \i + c a e\ 3 x) 
8c a g'i 2 xg^ 
1ic a g' x {Pi ~ P3) 
-2\{c a e\ 2 x + s w s a M w X 3 ) 

ejacScb 

1 c a eM T x 

2 M w s w ° ab 



1 eM T s a 

2 M\ySw 



1 e 



2 Cin Sf 



;7ac((l-2s 



2\0-r% 



9 „2 (l+7 5 )ct \ 



1 ieM T 5 

2 Mws w 'ab 



cb 



. 1 CgseVZ^ (1-7 5 ) 



2 s-u 



Toe" 



i/(l+7 5 )a6' 



- + 2s w s a3 iV%2/3- 

( Sa3 eM rX /2^ 



1 es a3 V2 „,/J (l-7 5 )cft 

2 ; 



-7 5 )ab 



2 Mn/s„ 
e^facScb 

2 M w s ro ° afe 

1 eMjjSg x 

2 M w s m 



2 Cin Si 



-7^((l-2s 



2s 



2 (1+7 5 W 



.1 jeM/±_ 5 

2 Mvi^Sii, *a& 



cb 



. 1 c ££ W2 ,M (1-7 5 ) 



2 



7ac" 



(c a2 eM / ,v / 2 (1 l 5)afe + 2s w s a2 M w y% 



2 M w s t 

1 es a2 V2 u (1^7 
"2 s... 7ac o 



(l+7 5 )aft ' 

2 
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Fields in 


the vertex 


Variational derivative ol Lagrangian by nelds 




vh\ 


w~ 


In/ (sa 2 eM u V2^t )ab 2s w c a2 M w y!< {1+ t )ab ) 


Sap 


Cbq 


W~ 


1 eV2V ca x -P (l-7 5 )ct 

2 s w "PI l ac 2 


Sap 


Cbq 


w~ 


1 ie ^ Vcs S vo (M s (1 -?f u M c {1+ t )ab ) 

2 M w s w "pq\^ rJ -s 2 c 2 ) 


Sap 


Sbq 


A, 


3 ^ u pq Iac u co 


Sap 


Sbq 


G 


ys A pq lab 


Sap 


Sbq 


Hi 


o \f Ms 5va$ab 


Sap 


Sbq 


n 2 


1 eM s s a x X , 

2 M w s w °pq°ab 


Sap 


Sbq 


z 


1 e x U (r> n 2 (l+7 5 )cb (o O „2 \ (l-7 5 )cb \ 


Sap 


Sbq 




1 ieM s x 5 

2 M w s w U P<1 lab 


Sap 


Sbq 


z' 


■^yi°pq JacVcb 


Sap 


t bq 


w~ 


1 ev^Vt, r (1-7 5 ) C 6 

2 s*, °P? ' ac 2 


Sap 


t bq 


UJ 


2 M w s ro °P?l iU s 2 iU * 2 J 


Sap 


Ubq 


w~ 


1 e^2V« s r M (l-7 5 )c6 

2 s w °P9 < ac 2 


Sap 


Ubq 


W 


lieM s V2V U sX (l"7 5 )a6 
2 M w s w U PQ 2 


tap 


hq 


w + 

A* 


1 eV2V ti)(i? u (1-Y°)cb 

2 s w U P1 < ac 2 


tap 


bbq 




lieV2V tb x ( t\j (l+7 5 )at (if (l-7 5 )at\ 
2 M w s M pg^ Vi & 2 iU * 2 J 


tap 


dbq 


w+ 

P 


1 eV2V td x M (l-7 5 )ct 

2 s w PI lac 2 


tap 


dbq 




lieAffVaVwi (1-7 5 U 
2 M ff s„ "P9 2 


tap 


Sbq 




1 ev^Vt, r (1-7 5 U 

2 U P<J '" c 2 


tap 


Sbq 


w 


2 AfwStu "PQV'^s 2 1 2 J 


tap 


tba 


A, 


^cupq jac u cb 


tap 


tbq 




9s\qlab 


tap 


tbq 


H x 


1 c a eM t x x 

2 M w s w °pq°ab 


tap 


tbq 


H 2 


1 eM t s a x x 

2 M w s w °pq°ab 


tap 


tbq 




4^W((3 - 44)^ - 44^) 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


tap 


tbq Z 


1 ieM t x _,5 

2 M w s w VI lab 




tap 


tbq Z /Lt 


-Ig'lSpqJacScb 




Uap 


ha W+ 


1 eV2V ub r u (1-7 5 ) C 6 

2 s w °pq < ac 2 




Uap 


hq W + 


lieM b V2V ub x (l+7 5 )ab 
2 M w s w °PQ 2 




Uap 


d bq W+ 


1 eV2V ud x „ju (1-7 5 ) C 6 

2 s w PI lac 2 




"Uap 




1 eV2V us x -P (l-7 5 )cb 

2 s w °P1 > ac 2 




Uap 


Sbq W + 


1 ieMsV^Vus x (l+7 5 )ai> 

2 M w s w U PQ 2 




Uap 


u bq 


3 ^pq^ac^cb 




Hap 


Ubq G ^ r 


ds^pqlab 




"Uap 


u bq %n 


-^«c((3-44)^- 


a„2 (l+7 5 )ct\ 
2 / 


Uap 


Ubq Z p 






w+ 


w~ z p 


-^[(pi-i$^-(p?-p£)<r 


-(p£-p£)s" P ] 


w+ 


w~ Z v 








w~ z 


4r-(^-^) 




C w+ 


c z w- 


A 4 




c w+ 


C z w~ 


— zeMvi/ 




c w+ 


C w ~ A, 






c w+ 


C w ~ H x 


1 c a eMw 
9 « 




c w+ 


C w ~ H 2 


1 eM w s a 

9 o 




c w+ 


c w - Z, 


c m e At 




c w+ 


C w ~ z 


1 ieM w 

2 S-u; 




c w+ 


c z w~ 






c w+ 


C z w- 


1 i(l-2a* )eM w 
2 




w + 


w- z v 






w + 


W~ z 
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Fields in the vertex 



Variational derivative of Lagrangian by fields 



Z C w b w V* * 

-epi 



C W ~ 

c w ~ 

C w- 
C w ~ 



w 



c z w+ 



<z 

c z 



c 



w+ 



c 



w+ 



c 



w+ 



Hi 
H 2 



C 



w+ 



Z, 



C 



w+ 



Pi) 



ieM\y 



1 c a eM w 
2 

1 eM w s a 

2 Sjn 



s«, Pi 

1 ieM w 

2 Sim 



C w- 


c z 


W + 


Cwe M 
s w Pi 


w- 


c z 


w + 


1 i(l-2s 2 w )eM w 

2 Cm Su; 


C z 


(jW+ 


W~ 


Cwe M 
s„ ^1 


C z 


(jW+ 


w~ 


1 ieMiy 

2 Cw S kj 


C z 


c w- 


W + 


Owe M 
Pi 


C z 


c w- 


w + 


1 ieM w 

2 CwSw 


C z 


C z 


H x 


1 c Q eA/w 


C z 


C z 


H 2 


1 eMvi/Sa 


C z> 


c z> 


Hi 




C z ' 


c z> 


H 2 


-2c Q M Zs _ i5 i 



1 Q 


vl\ 


w + 


1 c al eV2„p (l+7 5 )ct 

2 s w < ac 2 










vl\ 


w + 




+ 2s^s a iM H /y 1 ' 


(i- 


"7 5 )a6^ 


1 a 


vh\ 


W + 


1 es al y/2 u (l+7 5 )ci, 

2 s w " ac 2 








1 a 


vh\ 


w + 




- 2s u ,c a iM w y5 / 


(i- 


"7 5 )a6^ 


' a 


vl\ 


W + 


1 c a3 eV2 u (l+7 5 )ct 

2 s w > ac 2 








' a 


< 


w + 




+ 2s w Sa3M w y% 


(i- 




1 a 




W + 


1 es a3 V2 u (l+7 5 ) ci , 

2 s m " ac 2 
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Fields in the vertex 



Variational derivative of Lagrangian by fields 



A*S vl l 
Ma Vll 
Ma 



Ma 



Vl\ 



vht 



vll 



vll 



vll 



vll 



vll 
vll 

vll 



vll 
vll 
vll 
vll 
vll 
vll 



w 



w + 



w 1 

w + 

w + 
Hi 
H 2 



VL Zn 



vll 



vll 



vll 



vh\ 
vh\ 
vh\ 
vh\ 
vh\ 
vh\ 
vl\ 



vl 
vl 
vl u 
vll 



vll 



vh h 



vhl 



z' 
Hi 
H 2 

z 

z' 
Hi 
H 2 

z» 

z 

Z'» 
z' 

Hi 

H 2 



(s a3 eM T V2 {1+ l U - 2s w c a3 M w y 



2 M w s 

1 c a2 eV2 p (1+7 5 U 

2 s w < ac 2 



+ 2s w s a2 M w y v 2 



1 i_ 

2 M W i 

1 es a2 V2 p (l+7 5 )cb 

2 s,„ < ac 2 



(1-7 5 U' 



1 i_ 

2 M w a 

(caCaiVZxy^Sab + S a S a iM ul 5 ab ) 
^-{saCalV^Xy^ab ~ S a lC a M v 

2c w s w l°-cJ cb 

-ic a is a iV2y^ b 
-(1-2^)^7^ 
lib 



11 

" 2 x 



1 1 



((1 - 2s^)c Q V^^yi^aft + 2s Q S a iC a iM„ 

~((1 - 2s 2 al )s a \/2xy 1 (5 ab - 2s a ic a c a iM v 

1 c al es al a 5 

2 c„ S „ 7ac7 c fo 

-2c a is a i5 / 1 7 Q 1 c7cb 

x lab 

(c a c a2 V2xy^S ab + s a s a2 M V2 5 ab ) 
(s a c a2 V2xy!^5 ab - s a2 c a M U2 5 ab ) 

2c„s tl ,7ac7 c 6 
-ic a2 S a2 V / 2y2 7a6 

-(1-2^)^7^ 



'-Tab 



1 1 

"2 



1 1 

"2 x 



~((1 - 2s^ 2 )c a \/2a;y2^6 + 2s a s a2 c a2 M U2 5 ab ) 
((1 - 2s^ 2 )s a \/2xy^o6 - 2s a2 c a c a2 M v 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


< 




z. 


lc a2 es a2 u 5 
2 c w s w ' ac Icb 


< 


-hi 


z 


2 l \ 1 - ^ b a2)v ^i>2 lab 




vh\ 


z'» 


^a2 b a2y\ lac J c b 




vhl 


z' 


ic a 2M„ 2 s a 2 5 
x ^ab 




vl\ 


Hi 




< 


vl\ 


H 2 


So£(q r ^^/ 7 >'rii v ft L « M ft A 
^T\S a C a 3V AXV^Oab — S a ^C a lVl U3 a b) 


< 


vl\ 


z. 


2 c w s w < ac Icb 


< 


vl\ 


z 


<<<--ad € 'a.i V ^y% ! a b 


< 


vl\ 


z'. 


l x ^°aA>yi lac l c i) 


< 


vl\ 


z' 


iM "3 S a3„5 

x 'ab 


< 


vh\ 


H x 


— i-((f — 2s 2 n\c^\/~2xiin 5„h + 2s ni s„ic„iM,,„5„h) 

2 x V V do 1 a V i?3 &0 1 ^ a OL a ao^ao^^V^ au ) 




vh z b 


H 2 


2x \\ L z *a3/ 6 oV ^• L y3< J ab ^ b a3 L -a ( -a3 lvI us (J abJ 


ul 3 


vhl 


z» 


1 c aa es a[) u 5 

2 c w s w < ac Icb 


ul 3 
- l a 


vhl 


z 


2 l \ 1 - ^ b a3)v ^i>3 lab 


ul 3 


vhl 


z'» 


^ L -a3 b a3yi lac J c b 


< 


vhl 


z' 


iCasM^Saj, 5 




'V > 

x lab 


-h\ 


vh\ 


H x 


— £»i ( c , n \f 'Jt.i i 1 ? () „u — .s^r^i A/f,, ft„A 
rp y^ai^ct V ^• Xj y\ yJ ao "a u dl J - K - £ fi w aoJ 


-K 


vhl 


H 2 


_£al. ( a a -, ^F)'ri l v ft , -L r C i M ft A 

x \ b a b al\ ■^■ L y\'Jab ^ ( ^a.'^al lvl ui < Jabj 




vhl 


z, 


1 es li M 5 


-hi 


vhl 


z 


vy -al al V -^i/i 


-hi 


vhl 


z 'n 


V 1 ^ b a\)y\ lac l c b 


-K 


vhl 


z' 


a; 'ab 


-hi 


vhl 


Hx 


- £ f-(sa2C a V2xy%5 a b - S a C a 2M U2 8 a b) 


-hi 


vhl 


H 2 


- £ f-{s a Sa2V2xy%5ab + C a C a2 M U2 5 a b) 


-hi 


vh 2 b 


z. 


2c w s w 'acl ch 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 




-hi 


z 


ic~os~o \/2 / ii% '■y ■ 




-hi 


z'» 


^ x ■ 6d a2/yi l ac Icb 


-hi 


-hi 


z> 


ic2 a2 M "2 _,5 
x lab 


-K 


-hi 


H x 


x \°ao^ j av ^^y^^ab a ct t ~ j a6 ±vl v^ u ab) 


-K 


-hi 


H 2 


_£a3(„ „ ny/9Tii u /) i. 4- r r o M A 


-hi 


-hi 




2 c w s w ' ac Icb 


-hi 


-hi 


z 


^ao^flj v i a }) 


-hi 






(I - 2s 2 o)oW^ 5 , 


-hi 


-hi 


z' 


ic a3 M "3 -,5 
x 'ab 


K 


A v 


w + p w- a 


-e 2 (2q^ u qP (T - qWq ua - qP° ' q v P) 
^ \ y y y y y y / 




A v 


w + w~ 


2e 2 q pv 




Hi 


W+ w~ 


1 ic a e 2 uv 

2 s w y 




H x 


w+ W~ 


1 ic a e 2 uv 

2 Sw y 




H 2 


W+ w~ 


1 ie 2 s a uv 

2 Sw y 




H 2 


w+ W~ 


1 ie 2 s a uv 

2 Sw y 




W+ 


w- z a 


Cwe2 (2aV ja n v P n^ u aP a aPPa ua ) 
s w ^ y y y y y n j 




W+ 


w~ z 


2s w y 




w + 


W~ z 


2s w y 




w + 


w~ Z v 


(l-2s 2 w )e 2 U v 






G pr G as 


a 2 ( a^ p a ua L„+ - ' a vp L„+ + a p(J f™.* f„*+ 

ys\y y JpqtJrst y y JpqtJrst \ y y Jprtjqst 

_ pxj Q vp f f , a pv pa f f _ a fJ,p a va f f \ 
y y Jprtjqst ~ y y Jpstjqrt y y J pstj qrtj 


Ht 


Hx 


Hx Hx 


c 1 r A \ _i_ e 4 \ i „2 _2 \ \ 
— 0[C a Ax + S a A 2 + C a S a As) 


Hi 


Hx 


Hx Hx 


-3c a s a (2c 2 a Xx - 2s 2 a X 2 - (1 - 2s 2 a ) A 3 ) 


Hi 


Hx 


H 2 H 2 


-(644A1 + 644A 2 + (1 - 644)A 3 ) 


Ht 


Hx 


w+ w- 


1 4 e2 

2 q2 y 
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Fields in the vertex 


Variational derivative of La£ 


;rangian by fields 


Ht 


Ht 


w + w 


-(2c 2 \i + s 2 A^ 




Ht 


Ht 


Z u Z v 






Hi 


Ht 


z z 


-(2c 2 A! + s 2 A^ 




Ht 


Ht 


Z fj, Z' ' u 


8 s 2 a', 2 a^ y 




Ht 


Ht 


z' z' 


- (2s 2 \o + c 2 \*) 




Ht 


H 2 


H 2 H 2 


— 3c„ S^ ( 2s 2 Ai — 2c 2 Xo + (I 


- 2s 2 


Ht 


H 2 


w+ w~ 


1 c a e 2 s a Aiv 

z b w 




Ht 


H 2 


w + w~ 


— C S (2Ai — Aq"l 




Ht 


H 2 


Zfj, Z u 


1 c a e 2 s a u,v 

o r 2 „2 y 




Ht 


H 2 


z z 


— c^s^, (2Ai — Aq l 




Ht 


H 2 


Z n Z'u 


SCr^SrvQ^Q^ 




Ht 


H 2 


z' zf 






Ht 


w+ 


w~ Z v 


1 ic a e 2 uv 

2 c w " 




Ht 


w + 


w- z u 


1 ic Q e 2 uv 

2 c w y 




H 2 


H 2 


H 2 H 2 


_6f S 4 Ai + C 4 A? + C 2 S 2 \i) 




H 2 


H 2 


w+ w- 


1 e 2 s 2 




H 2 


H 2 


w + w~ 


-f2s 2 Ai + c 2 A^ 




H 2 


H 2 


Zfj, Z v 






H 2 


H 2 


z z 


-f2s 2 Ai +c 2 A s ) 




H 2 


H 2 


Z n Z 1 v 


8c 2 a', 2 ^ 




H 2 


H 2 


z> z> 


-f2c 2 A9 + s 2 A^ 




H 2 


W + 




1 ie 2 s a „uv 

1 r y 




H 2 


w + 


w- z v 


1 ie 2 s a uv 

2 c w y 




w+ 




w~ w- 


4 (2g^g pa - g^g u P - gfPg 




w+ 


w + 




o ,2 y 
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Fields in the vertex 



Variational derivative of Lagrangian by fields 



w+ w- z p z a 



W+ W~ z z 

W~ Zy Z 



W W w w 



w 



W~ Z u z 



w 1 w 



Z/j, Z v 



w'w z z 

w + w~ z' z' 

Zfj, Z v z z 

z z z z 



z z z' z' 



Z 11 Z' v z z' 



2c^y 

-4Ai 

2 Cw y 

1 (l-2slfe 2 gUU 



-2Ai 
-A 3 

1 e 2 ny 

-6Ai 
-A 3 
8g[ 2 g^ 
-6A 2 
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